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The purpose of this investigation is to demonstrate that 
the resistance of ceramic and ceramal materials to thermal 
shock can be determined by analyzing the nonsteady 
thermal stresses in the material. First the nonuniform 
temperature distribution is computed by using the heat 
conduction equation. This temperature distribution is 
then used to compute the thermal stress. The material 
parameters which enter into the calculations are the 
wefficients of thermal conduction and of thermal expan- 
sion, the ultimate strength, Young’s modulus, and Pois- 
on’s ratio of the material. These quantities together 
with the heating conditions then specify the resistance to 
thermal shock. The theory ig verified by comparing the 
results with NACA test data. 


gregh peter to thermal shock is the strength 
of the material against failure during rapid 
heating or cooling. The current interest in high 
lmperature designs using ceramic and ceramal mate- 
tals, which are relatively weak in thermal shock in 
mparison with metals, necessitates a closer exami- 
tation of this composite property. To help the search 
ior better materials, it is important to determine the 
physical and mechanical properties of the material 
which contribute to high thermal shock resistance. 
V. G. Lidman and A. R. Bobrowsky (1)? argued that 
ihe resistance to thermal shock is proportional to the 
talue of ook/aH, where k is the coefficient of heat con- 
luction of the material, a) the ultimate strength, a the 
vefficient of linear thermal expansion, and E, Young’s 
hodulus of elasticity. It is the purpose of this paper to 
ive a more complete theory of resistance to thermal 
tock, together with a comparison with experimental 
ata. 

The failure of material is the result of high stress. 
rapid heating or cooling, the stress is generated by 
he nonuniform temperature and hence nonuniform ex- 
insion. Then, the problem of resistance to thermal 
wek is really the problem of computing the thermal 
tess under such conditions. In order to correlate 
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the theory with experiment, the particular case ana- 
lyzed is that of a plate heated or cooled rapidly. For 
this simple shape, the computation is quite easy. In 
fact, the problem of thermal stress in a plate is already 
known (2), while the temperature distributions in the 
plate have also been calculated before (3). The spe- 
cific task of this discussion is then to bring these ele- 
ments together and properly interpret the results for a 
clear understanding of the problem of thermal shock. 


Thermal Stress in a Flat Plate 


Consider a flat plate with a uniform thickness much 
smaller than its lateral dimensions. Then if cooling 
and heating of the plate take place through heat trans- 
fer at the surfaces, the temperature of the material is 
different at different distances to the plate surface, but 
for the greater part of the plate the temperature is inde- 
pendent of the location in the plane of the plate. If 
the median plane of the plate is taken as the x, y plane 
(Fig. 1), then the temperature T is a function of z only. 
Let o,, and e,, be the normal stresses and 


ie b te b 


FIG. 1 DIAGRAM OF PLATE WITH BOTH SIDES UNIFORMLY HEATED, 
b HALF THICKNESS 
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normal strains in the x, y and z directions, respectively. 
The shearing stresses and shearing strains obviously 
vanish in this problem. The relations between the 
stresses and strains are 

= v( oy + o.)} + er 

ey = — oo, + os)} 
= Glee — + oy)} + er 


where e, is the thermal strain due to thermal expansion, 
and E and » are Young’s modulus and Poisson’s ratio. 
E and v may be functions of temperature. Due to the 
lack of shearing stresses and applied surface forces, 
a, is identically zero by the equilibrium of forces in the 
z direction. 

Now treat the case where the temperature is symmet- 
rical with respect to the median plane of the plate. 
Then the only deformation the plate can undergo is a 
uniform lateral expansion. Therefore, e, = e, and 
both are constants. Then o; = o,. Assume that the 
plate is not restrained at the edges, then the net tensile 


b 
force at the edges must vanish or, _, % dz = 0, 


where b is the half-thickness of the plate. By using 
these relations, it can be easily deduced that 


Loe 


b 
f 
—bv 
and hence 
b 
f d 
or = vy = 4 [3] 


These equations show that given temperature distribu- 
tion in the plate, the stress and strain at any point can 
be easily calculated, provided E, v, ey are known func- 
tions of temperature. If EZ, v are constants and take 
er = aT’, where T is the change in temperature from 
the initially uniform state and a the constant coefficient 
of thermal expansion, Equation [3] reduces to 


Ea {1 b 
= - 7) [4] 


By introducing the dimensionless parameters 


where 7 is any convenient reference temperature, 
Equation [4] simplifies to 


Hence the reduced stress is equal to the difference be- 
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tween the average reduced temperature and the local 
reduced temperature. 

For the more general case, where the temperature in J ¢ 
the plate is not symmetric with respect to the median 


plane but still independent of x and y, one need only I 
make the change ( 
a 

where e and 8 are two constants so far undetermined, ‘ 
0 


while the conditions ¢, = a, and o, = O remain true. 
That a linear variation of e, across the thickness is cor- 
rect follows from the symmetric consideration imposed 
by the assumptions on the uniformity of temperature Jy} 
with respect to x and y, and the absence of externa! con- 
straints. In addition to conditions of vanishing aver- 
age tensile stress across the thickness of the plate, there 
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FIG 2 DIAGRAM OF PLATE WITH ONE SIDE UNIFORMLY HEATED, 

b THICKNESS 
is the further condition of vanishing bending momeli{ji.t,j}, 
due to the assumption of free edges of the plate. heated 
The stress-strain relations are now 8 the 
lesiar 
Vv 

lary he. 


In view of later development in this paper, it is conVelfhyag: we 
ient to take the thickness to be b instead of 2b, a8 Hhivity aj 
the previous case. (Fig. 2.) Then the force allhng jr I 


moment conditions become mation 
b/2 
f o;dz = o| 
—b/2 | 
b/2 
ozzdz = 0 
—b/2 he boun 


By integrating now the first of Equations [7] with 
spect to z, and, by means of Equation [8], the follow! 
equations for the two unknown e and 8 are obtained, 
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b/2 b/2 
a+ Ez de = Eer 


b/2 b/2 b/2 Eerz 
= 


Here again, e and and hence e;, a7, o,, and e, are 
determined when the temperature distribution is known, 
and when E£, »v, and e, are known functions of tempera- 
ture. If e, = aT’, and if E, v, and a@ are constant, 
the following equation for the reduced stress is readily 


obtained. 
= T* + 12¢ 


where, as before, 


b/2 
€ 


Since for this case the thickness of the plate is 6, 
-(b/2) < z < (b/2). When T is symmetric with re- 
T*é dé = 0, then Equation [10] re- 
duccs to Equation [5]. 

Actually, to compute the thermal stress in a plate one 
needs to know the temperature distribution. For this 
purpose, consider two special cases: As the first case, 
take a plate heated or cooled uniformly on both sides 
by a constant heat source. As the second case, con- 
sider a plate uniformly heated or cooled on one side, 
and insulated on the other. As far as temperature dis- 
tribution is concerned, only the second case needs to be 
computed, since at the median plane of the plate in the 
first case there is no heat transfer due to the symmetry 
of the temperature distribution and the median plane 
tan be considered as insulated. This is the reason that 
the thickness in the symmetric case is taken to be 2b 
and the thickness in the second case is taken to be b. 


1/9 
spect to 2, 


Transient Temperature Distribution 


Proceed now to the calculation of the temperature 
listribution in an infinite plate with one side uniformly 
teated or cooled and the other side insulated. Take b 
is the constant thickness, and x and y axis of a Car- 
sian co-ordinate system in the middle plane. For 
mall strain, i.e., when the shape and size of the plate 
tre only slightly changed by thermal strain, the ordi- 
lary heat-conduction equation may be used. If k is the 
vefficient of heat conduction, h the surface conduc- 
ivity at z = —b/2, c the specific heat per unit volume, 

d if k, h, and ¢ are constants, the heat conduction 
quation is 

oT 
oz? 


the boundary conditions are 


oT b 
k oe = h(T = To), atz ~9 

oT b 

ae 0, atz = 2 


where 7) = constant is now taken as the temperature 
of the heat source. Initially the plate is considered to 
be at uniform temperature 7 = 0. The above equa- 
tions may be written in dimensionless form if one intro- 
duces 


[11] 
kt k 
f =o n= 


where R is the so-called resistance ratio, and f the 


Fourier number. Thus 

or* } 

and 

(12] 
and 

or* 1 

OE 0, até = 2 | 

T* 0,atf = } 


The solution to the above equation is well known (3). 
It is written down here for reference: 


= cos m(1/2 — £) 
(1 + R + sin 


1 


T* =1— [13] 


= 1 COS % 2 sin m 
where 7; are roots of the transcendental equation 


For small values of f, Equations [13] may be expressed 
asymptotically as (4) 


) 


T* = ~R+Re™ erfe O( 1/49) 


2R2f'/2 


f /2 Tt dé = 


f/R? Vf —(1/4f) 
Re erfe (v2) +0(f ) 


where erf(u) is the complimentary error function, de- 
fined as 


erfe(u) = 1 — erf(u) = ks. e—s? ds 


Computation of Thermal Stresses 


For the plate with both sides heated it follows from 
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FIG. 3 MAXIMUM THERMAL STRESS IN A LARGE PLATE WITH BOTH 

SIDES UNIFORMLY HEATED. AT £ = 0, TENSILE STRESS IN HEAT- 

AND COMPRESSIVE STRESS IN COOLING. AT § = +1, TENSILE 
STRESS IN COOLING AND COMPRESSIVE STRESS IN HEATING 


Equations [5], [13!, and [15] that 


k=0 i+ R + SIN 
1 
o* = —~R+ erfe — erfe 


erfe + ... [17] 


where now —1 <£ <1. It is to be noted that, since 


a T* — 7'* 


T* being symmetric and monotonic to either side of the 
z axis, maximum stress must occur at § = +1 and — = 
0. In the case of heating, i.e., T > 0, oat & = +1 
will be a compression and that at & = 0, a tension. 
The converse is true during cooling. Consequently in 
the following numerical computation, stresses are com- 
puted only até = —1l,andé =0 

Maximum tensile and compressive stresses in the 
plate are computed according to Equations [16] and 
[17], for various values of R. The result is plotted 
in Figs. 3 and 4. In the following paragraphs the ap- 
proximate formulas for maximum stresses at § = +1 
and ¢ = 0 for very small and very large values of R are 
given. 

By Equations [16] and [17], até = +1, 


9/4) 19) 


When R = 0, o* obviously reaches its maximum at 
f = 0. Hence it is expected that for very small R, 
Equation [19] should be adequate for computational 


FIG. 4 MAXIMUM THERMAL STRESS IN A LARGE PLATE WITH ONE 

SIDE UNIFORMLY HEATED. AT £ = 0, TENSILE STRESS IN HEATING 

ING AND COMPRESSIVE STRESS IN COOLING. AT £ = —!/2, TENSILE 
STRESS IN COOLING AND COMPRESSIVE STRESS IN HEATING 


purposes. ‘To determine maximum o%, set 


Equation [19] then gives 
1 2f'/2 


Hence, given ~/f/R, one can compute R and hence 
determine f. The applicability of Equation [20] is of 
course limited to small values of f since it is an asymp- 
totic solution. For four-decimal point accuracy, f is 
limited to values below 0.1. Corresponding to this 
value of f, R ~ 0.32. It can be seen from the structure 
of Equation |20] that, as R decreases to zero, V/5/R 
must increase to infinity, although f itself decreases. Bison 
In fact, hold, 


erfc VS = + 


rou 


To this order of approximation 


f= 3R? 


It can be shown that when R < 0.005, the above formula 
is adequate for four-place accuracy in f. When f is 
known, it is but a simple matter to calculate o*. Ac- 
cordingly as R—> 0, 


8 Ceramal 
+ Nis R's: [ 
: Titanium 
On the other hand, when R > ©, f also increases Sige. jin 
R increases. Hence, one may use Equation [18], theBilircon (z 
first two terms of which give accurate answers up t0 


four significant figures, when f is larger than 0.10. Magnesi 
After a similar analysis as outlined previously, one 0b N% sirec 
tains for maximum o* até = +1, ealciu 
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where mo, m are the first two roots of Ry tan n = 1, and 


1 1 +R + R%q? Rm? — 1 


1+ + R%? 1 — Ryo? 


Now, for large values of R, 


1 
1 


n = R- — | 
= [23] 
R-* + O(R-4) | 


— 


1 
n=rt 


By making use of these results, one obtains 


[1 + sap log OR + + 


or roughly, 


[24] 


This last result may be used when R > 10, to yield an- 
swers of two significant figures. It can be shown that 
when R > 1, Equation [24] already gives sufficiently 
accurate answers up to two significant figures. 

At the point = 0, preliminary calculations show 
that when R = 0, f = 0.102, for maximum o*. Hence 
the Fourier series solution [16] may be used for small as 
well as large values of R, when numerical solutions are 
not carried beyond four significant figures. While 
Equations [23] correspond to large values of R, one 
must use, for small values of R, 


+ R) + O(R?) 


m = + R) + OCR?) 


Accordingly, for x*at = 0, the following relations 
hold, 


Tmax 


32 (x — 2)% 


Equation [26] may be used when R > 3 for two signifi- 
cant figures. Comparison between Equations [22], 
[24], [26], and [27] shows that the ratio of maximum 
stresses on the surface to that at the center is about 2 
to 1 when R +=, and is about 10 to 3 when R > 0. 

For the plate with one side heated and the other side 
insulated, the maximum stresses are computed for § = 
—'/,andt=0. The results are plotted on Figs. 3 and 
4, It should be noted that in the present case the maxi- 
mum stresses at § = —!/. or £ = 0 may not be the ab- 
solute maximum stresses in the whole plate for given 
R. However, this deviation is not considered to be im- 
portant and the computed maximum stresses are rep- 
resentative. 

For comparison with the previous case, the following 
result is given for omax* at & = —'/: when R is very 
small, 


T 


Due to the difference in the definition in the thickness 
of the plate, for the same plate, FR in this case is only 
one-half of the R in the previous case. Equation [28] 
then shows that the thermal stress at the surface is 
smaller when only one side of a plate is heated, as is 
to be expected, since bending of the plate constitutes a 
relief of the thermal stress. 


Application to the Evaluation of Thermal Shock 
Resistivity of Certain Ceramic and Ceramal 
Materials 


An experimental method (1 and 5), currently used 
by NACA for the evaluation of the thermal shock prop- 
erties of ceramic and ceramal materials, is conducted 
in the following way: A specimen 2 in. in diam and 


| = (1 + log 12R + in. in thickness is subjected to repeated cycles of 
)| heating and cooling until fracture occurs. Heating 
or, roughly is performed in a furnace of controlled temperature, 
en (26a) and cooling in an air stream at about room temperature. 
a aces se The specimen was first subjected to 25 temperature 
(21) Band cycles with a furnace temperature of 1800 F. If it 
ula TABLE 
is ooo Physical ties 
k, Tensile Thermal shock evaluation 
(60-1100 F), Btu/(hr) strength, — (cycles completed) — 
Material in./in. F (ft?) (F/in.) E, psi psi 1800 F 2000F 2200F 2400F 
(22) Ceramal (80% TiC, 20% CO) 5.0 x 10-6 247 5.5 X 107 33,000 (1800 F) 25 25 25 25 
11,000 (2200 F) 
. Titanium carbide (TiC) 4.1 X 10-¢ 240 _ 5.5 X 10° 17,200 (1800 F) 25 25 25 21 
S$ Beryllium oxide (BeO) 5.28 X 10-6 233 4.28 X 10’ 6,200 (1800 F) 25 3 
the Zircon (ZrSiO,) 2.24 X 10-6 11.6 2.4 X 10° 8,700 (1800 F) 1 
to 14.5 (1832 F) 
),10. Magnesium oxide (MgO) 6.94 17.7(2192F) 1.24 10’ 3, 100 (1800 F) 
és 40.8 (2012 F) 
41% zirconium dioxide +6% 4.95 X 10-6 6.4(2192F) 2.5 x 10’ 6,750 


calcium oxide (CaO) 
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survived this treatment, 25 cycles were repeated in suc- 
cession with furnace temperature at 2000 F, 2200 F, and 
2400 F, except when failure occurred. The results ob- 
tained at the NACA Lewis Flight Propulsion Labora- 
tory, are collected in Table I. 

Test results showed that all specimens failed by ten- 
sion during. cooling, with the exception of stabilized 
zirconium dioxide, which failed by shear. Hence, only 
tensile strengths are given in Table I. The tensile 
strength for stabilized zirconium given in the table is 
considered close to the shear strength. 

W. G. Lidman and A. R. Bobrowsky reported (1) 
that the increase of thermal shock strength was accom- 
panied by an increase in the ratio koo/aH, where a» is 
the fracture tensile strength of the material under con- 
sideration. This correlation was qualitative and, in so 
far as the parameter koo/aF is not dimensionless, has 
only limited applications. However, the previous re- 
sults on the calculation of maximum stress in a plate can 
be used to determine a more satisfactory correlation 
between the thermal shock test results and the physical 
parameters of the material. 

The deciding factor whether a brittle material fails 
or does not fail is clearly the ratio of the failure strength 
to the actual stress. If the ratio exceeds one, it will 
fail, and vice versa. Moreover, general fatigue tests 
show that when a specimen is repeatedly stressed 
through cycles, the number of cycles that the specimen 
can survive bears a statistical relation to this ratio. 
This relation, however, may vary from material to ma- 
terial and becomes less definite when the stress cycles 
are operated close to the failure strength. Applied to 
the phenomena of thermal stresses, these considerations 
should still be true, if all changes in physical and me- 
chanical properties due to change of temperature are 
properly taken care of. Hence, one still can take this 
ratio as the criterion for thermal shock evaluations, and 
write 


where o now stands for the maximum thermal stress 
actually reached in the specimen. When S is less than 
unity, one can predict with certainty that the specimen 
will fail. Scan be called the thermal shock resistivity. 
To find the maximum stress in the circulate plate 
tested by NACA, the results of the previous section for 
a plate heated on both sides can be applied with the 
following assumptions: (1) The plate is infinite in ex- 
tent, and constant in thickness; (2) the physical 
parameters E, v, a, k, c, and h are constant; (3) the 
plate is heated (or cooled) uniformly on both sides by 
the same heat source (or sink). The time duration is 
such that at the end of each temperature cycle, a uni- 
form temperature is re-established in the plate; and (4) 
small strain. Applied to the present problem, assump- 
tions 4 and 3 are justified by the very nature of the ex- 
periment. Assumption 2 is made on account of the 
fact that little information is available on the variation 
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of E, v, a, k, c, and h with respect to temperature. As- 
sumption 1 seems warranted in view of the uncertainty 
introduced by assumption 2 in addition to the fact that 
the test specimen was reasonably thin. 

The theoretical results in the previous section com- 
puted for a plate heated or cooled on both sides may be 
expressed as 


where o is the maximum stress at Z = +b. These 
points are chosen simply because at these points, maxi- 
mum stresses occur both during cooling and heating. 
By the definition of S, Equation [29], 


Hence S depends only on R and another dimensionless 


t (1 
arameter ——_, —_. 

in general linear, as shown in Fig. 3. However, as 


The dependence of S on R is not 


when FR becomes sufficiently large, say R 
> 10, it is seen that S is then proportional to R or k/ 
hb. Thus, when plates of the same thickness are tested 
in identical ways, one has 
1 — v)ook 

(1 — 
for large values of R. This is the relation used by Lid- 
man and Bobrowsky. As noted before, this relation is 
sufficiently accurate for two significant figure computa- 
tions when R > 10. 

To compare the test results given in Table 1 with 
the theoretical thermal shock resistivity, h is taken to be 
the same for all specimens, since they were tested under 
identical conditions. The value h = 50 Btu/F ft? hr 
gives a satisfactory correlation. The room tempera- 
ture and v are taken as 60 F and '/;3, respectively. 
In the following table, the thermal shock resistivity 
theoretically evaluated is compared with the test 
results. 

In Table 2, except for ceramal, data given in Table | 
only have been used, together with Fig. 3 and Equa- 
tion [31]. For ceramal, at 2000 F, oo is taken to be 
22,000 psi, which is the mean value of the two tensile 
strengths at 1800 F and 2200 F. Where two different 
values for k are given, two values for S are evaluated. 
The reason that the tensile strength at high temperature 
is used instead of the tensile strength at low tempera- 
ture becomes clear, if one recalls that none of the speci- 
mens failed by compression during the test, and that the 
maximum stress occurs always on the surface. 

The assertion, that whether the specimen will or will 
not survive any thermal cycle will depend on the value 
of S, is satisfactorily confirmed by the correlation be 
tween experimental results and computations based on 
the simplified theory. The values of S for titanium 
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TABLE 2 


Thermal Furnace temperature 
shock —-Room temperature 60 F-— 


Material evaluation 1800F 2000F 2200F 2400F 
Ceramal (80% Exp. 25 25 25 25 

TiC, 20% CO) S 5.5 3.3 1.5 or 
Titanium carbide Exp. 25 25 25 21 


8 3.4 3.0 2.7 2.5 
Beryllium oxide Exp. 25 
1.17 1.06 
Zircon Exp. 1 
0.45 
0.54 
Magnesium oxide Exp. 1/, 
s 0.14 
0.28 
94% zirconium Exp. 0 
dioxide + 6% Ss 0.10 


calcium oxide 


carbide at 2000 F, 2200 F, and 2400 F are probably too 
high, since in the calculation the ultimate tensile 
strength at 1800 F is used instead of its exact values at 
these high temperatures, owing to the lack of data. 
The same comment goes for S at 2000 F, calculated 
for beryllium oxide. 


Concluding Remarks 


On the strength of this correlation with experiments 
presented in the previous section, one may conclude 
that for the type of experiment on thermal shocks 
described, the question of immediate failure upon ap- 
plication of a temperature gradient can be satisfactorily 
answered by calculating the thermal shock resistivity 
§ according to the present simple theory. When S is 
maller than unity, one can be quite sure that the speci- 
nen will suffer immediate failure. 

The question as to how many cycles a specimen can 
undergo before failure eventually occurs, remains open. 
The NACA test procedure described above cannot be 
wed as the basis of such an analysis, since the thermal 
tyeles were stopped as soon as 25 cycles were completed 
without failure. Yet, from the practical point of view 
his question is of greater interest than that of immedi- 
ate failure. Ceramic and ceramal materials are used 
gnerally for machine parts that have only a limited 
ile. Hence, in the study of the relationship between 
ihe number of cycles N and thermal shock resistivity, 
he primary interest is in low or mediumly high values 
{N, in contrast to the usual fatigue tests on metals. 
Merhaps an empirical relation between S and N of the 
bllowing type, 


N =em(S-1), §>1 


tay be constructed, where m is a positive number to be 
ktermined experimentally. This sort of relationship, 
i course, can be used only in a statistical sense. A 
eat deal of deviation is to be expected when S is close 
bunity. 

It may also be pointed out that the thermal shock re- 
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sistance is by no means an intrinsic property of a ma- 
terial. Thermal shock resistance depends rather 
strongly on the manner in which heat is supplied and 
the form of the specimen, as seen from the influence of 
b and h on S, according to the simple theory described. 
This fact limits greatly the usefulness of thermal shock 
experiments carried out on special specimens and under 
special conditions, to no more than an evaluation of the 
comparative merits of various materials against ther- 
mal shock. Such a test cannot predict the performance 
of any particular design subjected to specified heating 
or cooling conditions. For any specific design, the 
thermal stress induced in the material by the nonuni- 
form and unsteady heating and cooling must be first 
determined. Then with the known strength of the ma- 
terial under repeated stress, the performance of the de- 
sign can be ascertained. This is the rational theory of 
thermal shock resistance. It is the primary purpose of 
this paper to demonstrate what fundamental physical 
parameters are essential to the construction of a rational 
theory by which the thermal shock phenomenon may 
be analyzed for any given material of any shape and 
under any service condition. 
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This paper presents a summary of the results of a series 
of tests conducted to determine the tensile strengths of 
various types of steels when subjected to rates of heating 
comparable to those encountered by the metal parts of 
solid propellant rocket motors. The results of these tests 
indicate that time at temperature has a very definite effect 
upon the high temperature strength characteristics of the 
types of steels investigated under this program. These 
data indicate that the high temperature strength proper- 
ties based on relatively high temperature long-exposure 
periods, approximately 15 minutes, are conservative for 
design purposes when the exposure times do not exceed 5 
seconds. 


HE design engineer is continually attempting to 
improve the efficiency of rocket motors and at the 
same time accomplish a reduction in cost. One of 
the basic requirements necessary for achievement of 


rocket motor. 


Strengths of Several Steels for Rocket Chambers 
Subjected to High Rates of Heating 


By R. L. NOLAND 


Aerojet Engineering Corporation, Azusa, Calif. 


Institute, Columbus, Ohio. In the Aerojet tests the 
high temperature strength characteristics of various 
steels were determined under short-duration heating 
conditions similar to those encountered by the metal 
parts in solid propellant booster motors. The typical 
temperature gradient through the chamber wall at the 
end of propellant burning was duplicated in the test 
specimen. In the Battelle Memorial Institute tests the 
test specimen was uniformly heated through its cross 
section. Emphasis in the two programs was placed on 
the determination of an approximate strength-tempera- 
ture relationship for a large number of materials. 


Aerojet Test Apparatus 


Fig. lisa schematic sketch of an apparatus which was 
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these goals is a good working knowledge of the various designed and fabricated at Aerojet for the purpose of diti 
materials employed in fabricating the inert parts of a evaluating the strength of various steels when subjected der 
Accordingly, in an effort to aid the de- to heating conditions which are encountered in solid § *PP 
sign engineer in accomplishing his objectives, a material propellant rocket engines. This test apparatus con- be \ 
Hea 
PRESSURE GENERATING CHAMBER It 
PROPELLANT GRAIN TEST CYLINDER conc 
enco 
chan 
THERMOCOUPLES jecte 
the p 
ient | 
to th 
lowec 
exist 
1 
wall ; 
PRESSURE TAPS Mens 
NOZZLE 
FIG. 1 SKETCH OF MATERIAL EVALUATION TESTING APPARATUS 
where 
evaluation program was conducted, and the results of sists of a gas generating chamber where a solid propel# « 
this program are presented in this paper. The data lant charge is burned; a test cylinder where the ma hk = 
presented are mainly for short-duration units and were terial being studied is stressed and heated at a rap, _ 
obtained from tests conducted at Aerojet Engineering rate; and a nozzle section through which the generate} 7 _ 
Corporation, Azusa, Calif., and the Battelle Memorial gases are allowed to escape. This design was selecte@ h, = 
after considerable thought because of the flexibility ! 
1 Principal design engineer, Solid Engine Department. offers in solid propellant rocket materials testing workg—— _ 
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By suitable selection of the propellant charge and annul 
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FIG. 2 THERMAL GRADLENT THROUGH THE CHAMBER WALL AT THE 
END OF PROPELLANT BURNING FOR TYPICAL AEROJET BOOSTER 
UNIT 
as lar gas passage dimensions for a test series, heating con- 
of ditions imposed upon the test sample (the small cylin- 
‘ed der in Fig. 1) can be held constant over a wide range of 
id applied stress to the material, or heating conditions can 
on- | De Widely varied at a constant stress. 
Heating Conditions for the Test Cylinder 
It was desired to subject the test cylinder to heating 
conditions which are approximately the same as those 
encountered in solid propellant rocket engines. The 
chamber wall of a solid propellant rocket motor is sub- 
jected to a very high rate of transmission of heat from 
the products of combustion, resulting in a thermal grad- 
ient through the wall at the cessation of burning similar 
to that shown in Fig. 2. The analytical procedure fol- 
lowed in determining the heat transfer conditions that 
exist in a solid propellant rocket engine is as follows (1). 
1 The temperature at any position in the chamber 
wall at a given time increment is a function of two di- 
mensionless parameters. They are: 
at/Ly? = Fourier modulus............... [1] 
heLw/Kw = Biot modulus............. 
opel thermal diffusivity of the chamber material; ft®hr~! 
mag we = thermal conductivity of the chamber material Btu 
ft—*hr-! (°F 
raps t = time, hr 
rave" L., = thickness of chamber wall, ft 


8 


--MEAN CHAMBER WALL 


he = heat transfer coefficient, Btu ft~*hr-! °F-! 
The Fourier modulus is simply a measure of 


* Numbers in parentheses refer to Biblography at end of paper. 
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elapsed time, and the Biot modulus is a measure of the 
rate of heat transfer, both being in dimensionless form. 
The temperature at a given position in the chamber 
wall may be written in the dimensionless form: 


Tg — Tz _ unconsumed temperature difference at position 
Tg — 7; 7 in the chamber wall 


Tc = temperature of the gas in °F 
7; = initial temperature of the chamber wall in °F 

3 The dependence of the unconsumed temperature 
difference upon the Fourier and Biot moduli is very com- 
plicated when expressed analytically; therefore, graphs 
are used to show the relationship among the three par- 
ameters. For each position in the chamber wall, a 
separate graph is required. These graphs are used to 
determine the space mean, inside, and outside tempera- 
tures of the chamber wall. 


Aerojet Test Procedure 


Upon completion of fabrication of the test cylinder, 
measurements of wall thickness, circumference, and 
hardness were made at three positions along the cham- 
ber. At each position four wall thickness measure- 
ments were made, each at an angle of 90 deg to each 
other. A special tool was designed and fabricated for 
the purpose of obtaining accurate circumferential meas- 
urements. By making circumferential measurements 
before and after each test, it is possible to determine the 
amount of permanent set taken by the test cylinder. 
It must be remembered that these measurements are 
based on a gage length of approximately 14.4 in. (the 
circumference of the chamber) and, as a result, the 
effects of localized “necking down” of the chamber wall 
are minimized. 

Temperature data were obtained during testing by 
attaching thermocouples to the outside surfaces of the 
test cylinders at the three positions along the chamber 
where wall thickness and circumferential measurements 
had been made. The thermocouples were attached to 


FIG. 3 MATERIAL EVALUATION TEST APPARATUS 
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the chamber by resistance welding, eliminating the 
need for a foreign material which may affect the re- 
corded temperature measurements. 

During operation, pressure measurements were taken 
at the entrance and exit ends of the test cylinder. A 
photograph of the test apparatus in the firing bay ready 
for test is shown in Fig. 3. 

In order to obtain reliable strength-temperature data, 
it is desirable that the maximum pressure generated in 
the test unit occur at the end of propellant burning 
when the temperature is at a maximum. In Fig. 4 it 
can be seen that the maximum pressure occurs at the 
end of propellant burning. The recorded pressure-time 
curve for an actual test with conditions similar to those 
used in the theoretical calculation is also shown in Fig. 
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/ ~ 
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FIG. 4 THEORETICAL AND ACTUAL PRESSURE TIME CURVES FOR 
THE MATERIAL EVALUATION TEST APPARATUS 


The strength-temperature relationship of the mater- 
ials evaluated was determined by maintaining constant 
heating conditions and increasing the pressure progres- 
sively until rupture or excessive deformation occurred. 
This procedure was repeated at each heating condition 
until the strength-versus-temperature relationship was 
determined. After the strength-versus-temperature 
curve for the AISI-A4130 material had been determined, 
it was possible, by comparing the cold strength of a new 
material with that of the A4130, to predict an approxi- 
mate relationship between strength and temperature for 
the new material. Asa result, it was possible to obtain 
chamber rupture, or substantial deformation, at each 
heating condition without conducting a large number 


of tests. 
To determine whether a change in the structure of 


Material Carbon Manganese Phosphorus Sulphur 
AISI A4130 0.31 0.55 0.03 0.01 
N.E. 8630 0.29 0.82 0.02 0.02 
AISI C1015 0.11 0.49 0.015 0.03 
AISI 3140 0.39 0.73 Neg. Neg. 
AISI 2330 0.34 0.78 0.02 0.02 
Modified 


0.88 0.02 0.03 


AISI 1035 


1.03 0.02 0.03 


J-55 


156 


TABLE 1 CHEMICAL COMPOSITION STEELS TESTED IN THE MATERIAL EVALUATION TESTS 


Silicon 


the various materials had occurred during testing, some 
of the chambers were sectioned and metallographically 
analyzed. 


Aerojet Test Results 


A total of 69 material evaluation tests were conducted 
with various types of steels. The materials tested and 
the number of tests conducted with each material are 
as follows (the chemical composition for each of the 
steels is presented in Table 1): 

1 AISI-A4130 steel heat-treated to an ultimate 
tensile strength of approximately 150,000 psi—19 tests. 

2 Refire tests of chambers previously tested in 
(1)—four tests. 

3 AISI-A4130 steel heat-treated to an ultimate 
strength of approximately 190,000 psi—six tests. 

4 AISI-A4130 steel in the normalized condition— 
five tests. 

5 NE 8630 steel heat-treated to an ultimate 
strength of approximately 150,000 psi—seven tests. 

6 Boron-treated AISI-1035 steel, heat-treated to 
an ultimate strength of approximately 170,000 psi— 
seven tests. 

7 J-55 steel (National Tube Company material 
designation) in the normalized condition—-six tests. 

8 AISI 3140 steel, heat-treated to an ultimate ten- 
sile strength of approximately 170,000 psi—four tests. 

9 AISI C1015 steel, as purchased in the cold- 
drawn condition—five tests. 

10 Same as (9) above except that a 0.75-in. wide 
band of 0.125-in. thick mild steel was intermittently 
welded around the center of each chamber—two tests. 

11 AISI 2330 steel heat-treated to an ultimate ten- 
sile strength of approximately 150,000 spi—four tests. 

A typical set of test data is presented in Table 2. 
Graphical representations for the A4130 steel, heat- 
treated to 150,000 and 190,000 psi, and the N.E. 8630 
steel test data, are shown in Figs. 5 through 7, respec- 
tively, in order to demonstrate typical test data. Fig. 
8 presents curves of the ultimate strength versus tem- 
perature for all of the steels investigated in this pro- 
gram. For the 150,000-A4130 steel it will be noted that 
a large number of ruptures were not obtained, but there 
were sufficient near failures to establish the ultimate 
strength versus temperature curve. For instance, for 
practical purposes a 2.77 per cent permanent defor- 
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Test No.* 22 16 27 
Chamber material AISI AISI AISI 
A4130 A4130 A4130 

Initial strengh (Rockwell), psi X 10-* 128 153 153 
Final strength (Rockwell), psi X 10-* 123 153 153 
Applied stress, psi X 10-% 57 108 98 
Space-mean temperature, °F 930 950 1020 
Permanent deformation, %° 0.0 0.31 2.77 
Maximum pressure, psi 2560 4600 4220 
Wall thickness, in. 0.100 0.096 0.096 
Inside diameter, in. 4.36 4.36 4.36 
Plug diameter, in. ... 4.00 4.00 
Rate of temperature rise, °F sec 980 2100 1950 
Outside wall temperature at maximum 

pressure, °F 720 660 550 
Average pressure, psia ... 2762 2813 
Chamber no. 6 16 16 


* Test Series 45SC. ° R denotes rupture. 

mation, as obtained in Test No. 27, may be considered 
as closely approximating the ultimate strength of the 
material. The curves present approximately the nomi- 
nal values of ultimate tensile strength versus tempera- 
ture and should be satisfactory for design purposes in 
which the designs are based on the yield strength rather 
than the ultimate strength. It is customary in the case 
of steels to assume that the yield strength is equal to 
0.8 of the ultimate strength. A good correlation be- 
tween the allowable and applied stresses is obtained in 
comparing these data with those obtained in actual 
applications. This is clearly shown in Fig. 5, where the 
stress-temperature point for a short-duration rocket 
motor has been plotted. This unit, which had an 
original ultimate tensile strength of 150,000 psi, was 
stressed to approximately 90,000 psi, had a mean 
chamber wall temperature of 1040 F, and was perma- 
nently deformed to 0.5 per cent. 


AISI-A4130 Steel, Heat-Treated to an Ultimate Tensile 
Strength of Approximately 150,000 Psi 


For comparative purposes, Fig. 9 shows the ultimate 
strength versus temperature curves for the 150,000- 
A4130 steel based on the following data: 

1 The material evaluation tests. 

2 The Battelle Memorial Institute tests (see page 
161). 

3 Steel mill tests, in which the test specimens were 
heated to the test temperature for approximately 15 
minutes prior to testing. The strength versus tem- 
perature relationship formerly used at Aerojet was 
based on these data. 

4 The recommended tensile strength versus tem- 
perature curve for the purpose of the design of rocket 
motors is included in Fig. 9. The strength values 
shown on the design curve are 80 per cent of the ulti- 
mate strength based on the material evaluation tests. 

In Fig. 9 it will be noted that a good correlation is ob- 
tained between the material evaluation data and the 
Battelle Memorial Institute data. As discussed on 
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TABLE 2 SUMMARY OF EXPERIMENTAL DATA, CHAMBER MATERIAL EVALUATION TESTS 


29 47 36 43 44 50 51 53 56 
AISI NE NE NE NE NE NE NE 
A4130 8630 8630 8630 8630 8630 8630 8630 1035 
153 141 142 141 144 145 147 143 171 
... ‘141 115 140 144 96 
137 30 78 74 38 86 84 50 55 
600 1300 1200 1140 1240 880 1000 1270 1070 
R R R 0.69 0.09 0.05 0.0 R 0.30 
5920 2000 3300 3320 1660 3730 3750 2160 2630 
0.096 0.094 0.093 0.098 0.095 0.098 0.098 0.094 0.091 
4.36 4.36 4.36 4.386 4.36 4.36 4.36 4.36 4.18 
4.0 3.32 3.32 4.0 3.32 
900 1270 1300 1100 1400 820 1070 =1210 =1100 
440 1125 1080 980 1050 750 805 1115 940 
2517 1442 2791 2920 2140 
16 23 39 25 24 31 32 35 56 


page 161 of this paper, the specimen is uniformly heated 
through its cross section in the Battelle Memorial Insti- 
tute tests, whereas in the material evaluation tests there 
is a temperature gradient through the chamber wall at 
the end of propellant burning. These data indicate 
that for AISI-A4130 steel, quenched and tempered to 
an ultimate tensile strength of 150,000 psi, the tempera- 
ture gradient through the chamber wall may be used as 
a criterion for determining the strength of the chamber 
at the end of propellant burning. 
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ULTIMATE TENSILE STRENGTH | 
BASED ON THE BATTELLE 
MEMORIAL INSTITUTE DATA 


ULTIMATE TENSILE 
STRENGTH BASED ON 
EVALUATION TESTS 


(20+ RECOMMENDED TENSILE 


STRENGTH DESIGN ie \ \ 
\ 
\ 
\ 
\ 


| 


2 | DESIGN CURVE PREVIOUSLY \ 

i | USED FOR DESIGN PURPOSE, 1 \ 

100 \ 

4 \ 

\ 

N \ 
\ 

80 \ 
ULTIMATE TENSILE STRENGTH \ 
PREVIOUSLY USED FOR DESIGN \ 
PURPOSES 

\ 
\ 
\ 
\ 
50 \ 
200 400 1000 1200 + 1400 


600 800 
TEMPERATURE, °F 


FIG. 9 TEMPERATURE VERSUS TENSILE STRENGTH FOR AIsI-A4130 
STEEL HEAT-TREATED TO ULTIMATE STRENGTH OF 150,000 PsI 
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A comparison of the material evaluation data and 
that previously used for design purposes (see Fig. 9) 
indicates that time at temperature has an appreciable 
effect upon the high-temperature-strength properties of 
the steel. In the material evaluation tests the speci- 
mens were heated for less than two seconds, whereas 
the other data are based on “long time”’ heating of the 
test specimens. The difference in strength properties 
obtained by the two methods of test is probably due to 
the extent to which the grain structure of the steel is 
altered prior to, or during, testing. For the “long per- 
iod’”’ tests a structural change probably occurs prior 
to test, whereas the material evaluation tests were con- 
ducted in such short periods of time that it was not pos- 
sible for a grain structure change to occur. A 150,000- 
A4130 chamber was quenched with water approximately 
five sec after completion of propellant burning, and a 
metallographic examination of the steel following the 
test indicated that no structural change had occurred 
even though it was at a temperature of 1180 F for at 
least five sec prior to quenching. However, on a basis 
of the other tests, it is known that if the mean tem- 
perature exceeds approximately 950 F (the original 
tempering temperature of the chamber) for an appre- 
ciable time and the chamber is allowed to cool slowly, 
a structural change of the steel occurs. This indicates 
that time at temperature is an important factor in alter- 
ing the structure of the steel. If the mean temperature 
of the chamber wall remains below the original temper- 
ing temperature, the strength of the chamber wall is 
not affected, with the exception of a slight amount of 
alteration that might occur at the inner surface of the 
chamber. 


AISI-A4130 Steel, Heat-Treated to an Ultimate Tensile 
Strength of Approximately 190,000 Psi 


In Fig. 8 it will be noted that at temperatures above 
1000 F there is little difference between the strengths 
of the A4130 steel heat-treated to 150,000 psi and 190,- 
000 psi. As a result, it may be concluded that there is 
little advantage in using the higher strength steel in 
applications when the mean chamber wall temperature 
exceeds 1000 F. The greatest advantage of the 190,- 
000-A4130 steel? is to be gained by maintaining its 
mean chamber wall temperature at values below 1000 F. 
It may be advantageous to use the higher strength steel 
if the metal part is to be subjected to repeated firings, 
because after test it retains a higher ultimate strength 
than the lower strength steel. 


AISI-A4130 Steel in the Normalized Condition 


The data in Fig. 8 indicate that the strength of the 
44130 steel in the normalized condition does not de- 

§ 190,000-A4130 will henceforth be used as an abbreviated form 
of referring to AISI-A4130 steel heat-treated to an ultimate ten- 
tile strength of approximately 190,000 psi. 
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crease as rapidly with increasing temperatures as the 
quenched and tempered steels. Because the quenched 
and tempered steels lose strength so rapidly at the high 
temperatures, the data tend to indicate that the nor- 
malized steel will function as satisfactorily as the 
quenched and tempered steels at temperatures above 
1400 F. 


Cold-Drawn AISI 1015 Steel 


Fig. 8 indicates that C1015 steel is good only for 
applications at high temperatures where low magnitudes 
of stress are involved. However, if the temperature of 
the steel is maintained below 600 F, it can be used 


effectively. 


J-55 Steel 


The J-55 steel in the normalized condition has tem- 
perature-strength properties very similar to those of 
the C1015 steel (see Fig. 8). One possible advantage of 
this steel over the C1015 steel may be that it can be 
welded continuously and then renormalized and have its 
original strength properties restored. 


N.E. 8630, Heat-Treated to an Ultimate Tensile Strength 
of Approximately 150,000 Psi 


The results of the tests conducted with this steel indi- 
cate that its strength-temperature relationship is al- 
most identical to that of the 150,000-A4130 steel. 


AISI 2330 Steel, Heat-Treated to an Ultimate Tensile 
Strength of Approximately 150,000 Psi 


The results of the material evaluation tests conducted 
with this material indicate that it has strength proper- 
ties equal to approximately 95 per cent of those for the 
150,000-A4130 steel at temperatures above 1000 F. 


AISI 3140 Steel, Heat-Treated to an Ultimate Tensile 
Strength of Approximately 170,000 Psi 


The results of the material evaluation tests con- 
ducted with this material indicate that in the tempera- 
ture range of 1100-1200 F the AISI 3140 material has 
strength properties equal to approximately 88 per cent 
of those of the 150,000-A4130 steel. 


Boron-Treated AISI 1035 Steel, Heat-Treated to an Ulti- 
mate Strength of Approximately 170,000 Psi 


This material was originally purchased as AISI 1035 
steel, and chemical and physical tests conducted by 
Aerojet indicated that it conformed to the specifications 
for AISI 1035 steel. However, when this steel was oil- 
quenched and then tempered at a temperature of 1000 
F, it had a Rockwell “C” hardness of approximately 37. 
Because normal AISI 1035 steel has a Rockwell “C” 
hardness of approximately 30 following an oil quench (2), 
a more thorough chemical analysis of the steel was con- 
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ducted to determine whether sufficient residual consti- 
tuents were present to give the steel its high harden- 
ability qualities. The residual constituents, such as 
nickel, molybdenum, chromium, and copper found in 
the steel, are not in sufficient quantity to give the steel 
its high hardenability characteristics. However, as 
presented in the Transactions of the American Society 
for Metals, Vol. 37, 1946, it is known that very small 
quantities of boron (approximately 0.0025 per cent) 
have an appreciable effect upon the hardenability char- 
acteristics of a plain carbon steel (3). Asa result, a spec- 
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HEATING RATE INDICATED ON CURVES 


trographic analysis of the steel was conducted, and it 
was found that it contained approximately 0.002 per 
cent boron. Consequently, it was possible to include a 
boron-treated steel in this evaluating program. The 
results of this study are indicative of the potential use 
of plain carbon steels treated with boron, because the 
hardenability characteristics of the steel studied in this 
program were slightly superior to those of AISI-A4130 
steel. It is realized that this may not be a true boron 
steel, where the boron is added to the metal under care- 
fully controlled processing techniques. 
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This boron-treated AISI 1035 steel is approximately 
85 per cent as strong as the 150,000-A4130 steel at tem- 
peratures above 1000 F. The boron-treated steel ap- 
pears to lose strength with increase in temperature 
more rapidly than the A4130 steel. Because of the high 
strength characteristics of this material at temperatures 
below 900 F, it would be advantageous to utilize this 
steel in applications in which the temperature did not 
exceed 900 F. 


Battelle Memorial Institute Test Results 


Tensile tests with high rate of rise in temperature 
were conducted by Battelle Memorial Institute in order 
to obtain a comparison between specimens uniformly 
heated and those with temperature gradients through 
their cross sections. In addition to obtaining tempera- 
ture-strength data by this method of testing, it is pos- 
sible to obtain the variation of modulus of elasticity 
with changes in temperature and stress. These data 
are invaluable for use in thermal stress problems. 

A brief description of the test setup employed by 
Battelle follows: 

| The test section of the specimen is '/2 in. wide and 
approximately 8 in. long. The over-all length of each 
specimen is approximately 17 in., and the width at the 
ends where the specimen is clamped during test is 2 
inches. 

2 Heating is obtained by passing a high current 
through the specimen from an electric are welding trans- 
former. It is possible to vary the rate of rise in tem- 
perature from 10 F/sec to 1300 F/see. 

3 Two test specimens clamped together are tested 
simultaneously. The rise in temperature is measured 
with thermocouples clamped between the two speci- 
mens. 

4 The tensile load is applied to the specimen by 
means of a dead load through a 9:1-ratio lever arm. 
The test equipment does not allow for a rate of change 
in the applied load, and, as a result, it is necessary to 
conduct the tests under constant loading conditions. _ 

Tests were conducted with Uniloy 19-9DL stainless 
steel and AISI-A4130 steel, quenched and tempered to 
an ultimate strength of approximately 150,000 psi. 


Uniloy 19-9DL Stainless Sieel 


The results of tests conducted with the 19-9DL stain- 
less steel specimens in a 1200 F stress-relieved condition 
are shown in Figs. 10 and 11. The figures are repro- 
ductions of the curves forwarded to Aerojet from the 
Battelle Memorial Institute. The termination of each 
curve indicates where failure of the specimen occurred. 
In Fig. 11 it will be noted that a variation in the rate of 
temperature rise from 280 F/sec to 550 F/sec has little 
effect upon the tensile strength properties of the ma- 
terial. 

By use of the data from Fig. 10 and assuming that 
there is no strain at zero stress, it is possible to plot con- 
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stant temperature curves of applied stress versus per 
cent deformation as shown in Fig. 12. These curves 
represent a conservative average of the plotted points 
and are considered to be only approximate because of 
the lack of sufficient data. However, if it is assumed 
that these are valid stress-strain curves, it is possible to 
ascertain values of tensile strength versus temperature. 
These data are graphically represented in Fig. 13. 
Fig. 13 also includes curves of the yield (0.2 per cent 
permanent-set method) and ultimate tensile strengths 
based on short-time mill test data for comparative pur- 
poses. These data indicate that for very short periods 
of heating the 19-9DL stainless steel has considerably 
higher tensile strength than is shown by the short-time 
mill test data. 

By comparing these data with those obtained for the 
A4130 steel, as shown in Fig. 9, it will be observed that 
the 19-9DL stainless steel, in the 1200 F stress-relieved 
condition, appears to be superior to the AISI-A4130 
steel, quenched and tempered to an ultimate tensile 
strength of 150,000 psi. However, with the A4130 
steel, it is possible to obtain a weld strength which is 
equal to, or greater than, the strength of the parent 
metal. When the 19-9DL stainless steel is welded, the 
areas adjacent to the weld lose considerable strength. 
Accordingly, in many applications the A4130 steel will 
actually be superior to the 19-9DL stainless steel. 


AISI-A4130 Steel, Heat-Treated to an Ultimate Tensile 
Strength of 150,000 Psi 


The graphical representation of the data on defor- 
mation versus temperature received from Battelle 
Memorial Institute is shown in Fig. 14. The ultimate 
strength versus temperature based on these data is 
shown in Fig. 9. By comparing these data with the 
data obtained with the test cylinders, it will be observed 
that the temperature gradient through the cross section 
of the specimen does not have a great effect upon the 
strength properties of the material. This is for AISI- 
A4130 steel with an approximate ultimate tensile 
strength of 150,000 psi. The temperature gradient 
may have an appreciable effect upon the strength of 
other materials, for it is a function of the relationship 
among stress, strain, and temperature. 
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Aspects of Combustion Stability in Liquid 


Propellant Rocket Motors 


Part I: Fundamentals. 


Low Frequency Instability With 


Monopropellants 


By L. CROCCO! 


Daniel and Florence Guggenhe'm Jet Propulsion Center, Princeton University, Princeton, N. J. 


Existing theories are able to explain part of the observed 
effects of low-frequency unstable combustion with the 
assumption of a constant time lag between the propellants 
injection and their transformation into hot gases. A dis- 
cussion shows that the time lag cannot be constant if the 
pressure oscillates, but must be correlated to the pressure 
history during the time lag itself. A reasonable formula is 
suggested to represent this correlation, and this formula 
is applied to the determination of the critical cond ‘tions. 
It is first shown that instability can ar‘se even for con- 
stant rate of injection. Secondly, it is shown how the 
variation of the injection rate affects the critical conditions 
for monopropellant feeding systems of constant rate or 
constant pressure type, having given characteristics of 
inertia and elasticity. Finally, the influence of non- 
uniform values of the time lag for different portions of the 
injected propellant is discussed. The bipropellant case 
and the case of high frequency instability will be discussed 
in Part II, to be published in next issue of the *‘ Journal of 
the American Rocket Society.” 


Nomenclature 


Combustion Chamber 


{ = time 

instantaneous value of the time lag 

value of the time lag in steady operation 

= pressure exponents of pressure dependence of the 
processes taking place during the time lag 

instantaneous pressure in the combustion cham- 
ber 

pressure in the combustion chamber in steady 
operation 

(p — p)/p = fractional variation of pressure in 
the combustion chamber 

instantaneous rate of injection, burning and ejec- 
tion of propellants 

= common value of these in steady operation 
= (m; — m)/m; pw, = (m, — m)/m = fractional 
variation of injection and burning rate 
A = instantaneous mass of gases in the combustion 
chamber 

F = same in steady operation 

f = M,/m = gas residence time in steady operation 

6, + + = total residence time of propellants in 
steady operation 

characteristic velocity and length of the rocket 


T = 


t = 


tity, m, 


, L* 
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R, y 
Ty Pg 


z= t/0,;6 


> 


Subscript 


gas constant and adiabatic index of combustion 
gases 

absolute temperature and density of combustion 
gases 

7/0, = reduced time and reduced time lag 

\ + iw = root of the characteristic equation with 
the reduced time as the independent variable 

reduced amplification coefficient 

reduced angular frequency 

angular frequency 

period of oscillations 


critical conditions of incipient instability 


Monopropellant Feeding System 


Pr 


= (pi — pi)/2Ap = 


instantaneous pressure at that place in the feed- 
ing line where the capacitance representing the 
elasticity is located. 

same in steady operation 

p.1 — p = injector pressure drop in steady opera- 
tion 

relative variation of p, 

regulated gas pressure for constant pressure 
supply 

area of the injector ports and of the cross section 
of the feeding line 

velocity of the fluid in A; and A; 

density of the propellant 

variable volume of the capacitance introduced in 
the feeding system to represent its elasticity 

dC/dp = compressibility coefficient 

mass of propellant in the line 

length of the line 

fractional length for the constant pressure supply 

2px Ap/ mo, = elasticity parameter of the line 

Im/2ApA4, = inertia parameter of the line 

p/2Ap = pressure drop parameter 

instantaneous mass flow and its fractional varia- 
tion upstream of the capacitance 


Nonuniform Time Lag 


f(0 < f < 1) = location of a given fraction of propellant in the 


a(f) 
5 


df) 


jet from the injector 
time lag corresponding to the fraction considered 
corresponding reduced time lag 
average reduced time lag 
5(f) — 5m = variation from the average 
extreme values of ¢(f) 
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1 Introduction 


HE pressure in the combustion chamber of a 

rocket motor is determined by the balance be- 
tween the mass flow of propellants into the chamber 
from the feeding system and the mass flow of gases 
ejected from the combustion chamber through the 
nozzle. If the feeding system has the desired properties 
of stability, that is, is able to provide under constant 
operating conditions a constant flow of propellants 
through the injector ports, then, since the exhaust 
through the nozzle is also constant for constant opera- 
ting conditions, the pressure in the combustion cham- 
ber would have to stay perfectly constant as long as no 
voluntary change is introduced in the operating par- 
ameters. Direct measurement of the pressure shows 
that this constant pressure is never exactly obtained, 
and that even when the combustion looks smooth, 
small, high-frequency fluctuations of the pressure are 
always present. But what is more interesting for the 
practical engineer is the case of rough combustion, 
when the pressure fluctuations are of large amplitude 
and can produce unfavorable effects on the operation 
of the rocket and eventually on its life. 

It is only in recent years that the observations on 
rough combustion have been co-ordinated and a dis- 
tinction made between low-frequency oscillations 
(chugging: generally below 100 cycles per second) 
and high-frequency oscillations (screaming: many 
hundreds of cycles per second). Also, a theory has 
been advanced and studied quantitatively (1, 2, 3)? 
which appears to explain satisfactorily some of the ob- 
served phenomena of the low-frequency instability. 

The purpose of this paper is to extend the discussion 
of the existing theory, and to introduce the considera- 
tion of a new factor of instability which not only may 
affect considerably the results in the low-frequency 
case, but which affords an explanation of the high- 
frequency instability. 


2 The Time Lag: Physical Facts 


The key to the existing theory resides in the con- 
sideration that the processes which are responsible for 
the transformation of the propellants into high tem- 
perature gases require a certain time. That the exist- 
ence of this time lag is a possible reason for instability 
can be qualitatively seen as follows: If the development 
of high-temperature gases follows immediately the in- 
jection of the propellants, an accidental increase in the 
chamber pressure would result in, on one side, an in- 
crease in gas flow from the rocket, and on the other 
side, a decrease in propellant flow to the rocket chamber 
(if the feeding system is sensitive to pressure vari- 
ations) and an immediate decrease in burning rate, 
so that the accidental over-pressure would drop to zero 
with time, and a self-excited oscillation would not be 


2 Numbers in parentheses refer to the Bibliography at end of 
Part I. 
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possible. But if, on the contrary, the burning does not 
immediately follow the injection, then a self-excited 
oscillation can exist because a reduction in mass flow 
due to an excess in chamber pressure results in a de- 
creased burning rate at a later time, when, provided 


the time lag and the period of the oscillations are in a 
certain relation, there is a pressure defect. 

The physical and chemical processes taking place 
during the time lag are quite complicated and, in any 
case, are not sufficiently known, so that a quantitative a 
treatment is impossible today. Nevertheless, some p 
qualitative information can be obtained from the fol- la 
lowing considerations: 

(a) Let us consider first the simpler case of a mono- al 
propellant. As soon as the liquid leaves the injector ec 
in the form of one or many jets, the combined actions m 
of the relative velocity with respect to the gases and al: 
of the capillarity break the jets into droplets. The no 
rapidity and the fineness of the resulting atomization siz 
depend substantially on the type of injector and on the cal 
density of the gases, and generally increase with the tal 
latter, so that for a given injector an increase in cham- J pel 
ber pressure has the effect of producing faster and finer J the 
atomization. Once the droplets are formed they move §f vay 
through a mass of hot gases, produced downstream by § cre: 
the previously burned propellant and recirculated § the 
backwards by eddies. This recirculation is always pres- §} con 
ent in a rocket chamber and is believed to be an essen- § con: 
tial part of the mechanism of combustion. Its intensity §f bip: 
will depend substantially on the type and distribution § the 
of the injectors, but, except for eventual periodical § lag 
variations in the flow pattern, can be considered to be § pres 
constant for a given injector system. The relative § prop 
motion of the droplets with respect to the high-tem- § in lic 
perature gases will produce heating and partial vapori- § tion 
zation of the droplets, and the rapidity of this process 
increases with the temperature difference (which is 3 7 
fixed) and with the density of the gases. Also, it will in- Fr 
crease with decreasing size of the droplets. For both Jfa sm 
reasons we see that an increase in pressure produces 4 {small 
faster heating of the droplets. So far, the process is § mixec 
only physical. The rest, which involves chemical ffnal 
transformations, can take place in two different ways. §tinuit 
Following the first in the heated body of the droplet Bitically 
(and eventually of the formed vapor) some exothermic fof the 
reactions begin to develop at a rate that increases with fivery s 
the temperature, so that once they have started they fiore, t) 
are also accelerated, and in a very short time a com- fitime a 
plete transformation is obtained into the final high- ffinto he 
temperature products. The rapidity of these reactions iific y¢ 
will also in some cases (especially when they take place Bihat of 
in the vapor phase) increase with the pressure. The ihe enc 
second way in which the transformation can happen is fhellant: 
by a surface process, analogous to the deflagration of flicket 
solid propellants, on the surface of the droplets.’ Here § A se, 

’ The evidence of such possibility for mixed bipropellants is lime lag 
afforded by recent experiments of Tait, Whittaker, and Williams xactly 
(4). llant, 
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again, as the rate of deflagration increases and the size 
of the droplets decreases with the pressure, the time 
necessary for the transformation in hot gases will de- 
crease with the pressure. Finally, we can conclude 
that an increase in chamber pressure produces a de- 
crease ‘n the partia' times required for the different 
processes to take place, and therefore a decrease of 
their sum, which represents the time lag. 

(b) The bipropellant case is more involved because 
a necessary prerequisite for any combustion to take 
place is a certain degree of mixing of the two propel- 
lants. In this case, no matter how finely and quickly 
each one of the propellants has been atomized, heated, 
and vaporized, nothing can happen until it comes in 


r contact with the other propellant. Now the rate of 
8 mixing depends essentially on the injection system, and 
d also to a smaller extent on recirculation, but probably 
e not on pressure. However, the conditions, droplet 
n size, temperature, vaporization, and eventually chemi- 
e cal transformation, of the propellants when the mixing 
ie ff takes place depend considerably, as in the monopro- 
- pellant case, on the chamber pressure; and therefore 
ar § the rapidity of the following vapor-vapor reaction or 
re § vapor-liquid reaction will also be increased by an in- 
yy § crease in pressure not only because of a direct effect on 
“1 @ the kinetics of the reactions, but also because of the 
s- § conditioning of the propellants before mixing. The 
n- § consequence of these considerations is that also in the 
ty § bipropellant case the total time lag will decrease when 
on § the pressure increases, though that part of the time 
‘al @ lag that represents the mixing delay is unaffected by 
be § pressure. (An exception is afforded by two self-igniting 
ve § propellants mixed in a nonatomized form and reacting 
m- § in liquid-liquid phase, a case in which a pressure varia- 
ri- § tion will have but negligible effect on the time lag.) 

yr 3 The Time Lag: Approximate Assumptions 
in- From what has been said, it follows that if we track 
oth Ja small amount of monopropellant, or corresponding 
s& §small amounts of the two propellants destined to be 
; is §mixed and burned together, from the injection to the 
cal Bfinal conversion in hot gases, though no real discon- 
8. tinuity in the physicochemical processes exists, prac- 
let F tically all of the chemical reactions, and the largest part 
mic fof the temperature and volume increase, take place in a 
ith Bvery short time near the end of the time lag. There- 
hey fore, the time lag can be approximately considered as the 
om- fiime after which the liquids are suddenly transformed 
igh- Binto hot gases at the final temperature, and as the spe- 
jons Bcific volume of the liquids is negligible with respect to 
lace Hihat of the gases (even at large pressures), it is only at 
The fthe end of the time lag that the introduction of the pro- 
nis fhellants will be felt in the balance of masses in the 
n of fhocket chamber. 

{ere # A second approximation consists in taking the same 
ts is ime lag for all the particles of propellants. This is not 
jams @actly true even in the simple case of the monopro- 


llant, since the conditions of the particles are not uni- 
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form in all the spray; and it is certainly a very rough 
assumption in the bipropellant case since the mixing 
delay is fundamentally different for different particles 
injected simultaneously. In Section 9 we shall discuss 
the effect of the nonuniformity of the time lag. For the 
moment, however, let us accept the simple assumption 
of a uniform time lag, that is, all the partic'es injected 
simultaneously are burned simultaneously. 

The actual determination of the time lag for a given 
combination of propellants, given designs of injecting 
system and rocket chamber, and given pressure level 
in the chamber, can be reached only experimentally. 
If, without changing the other conditions we alter the 
combustion pressure, the time lag, as it has been seen, 
must change, and the variation should also be deter- 
mined experimentally. As such experimental informa- 
tion is not yet available, it is necessary to assume some 
plausible and simple law of variation. For instance, 
one could assume that the time lag is inversely pro- 
portional to the pressure or, more generally, to some 
power of the pressure, so that 


Tp" = 7p" = const 


where the superposed bar indicates some fixed refer- 
ence condition. A more general formula could be 
tn; 

T = 7) = const; mp" = 7p” = const 
For small pressure variations it is possible to approxi- 
mate both formulas with a linear relation, and the two 
become equivalent, with certain relations between the 
constants involved. 

So far, the combustion pressure has been supposed 
to be constant during the operation, and only the effect 
of changing this pressure level has been considered. 
What if the pressure is oscillating? In this case the 
pressure will change during the time lag, and it is its 
integrated effect during the time lag which determines 
the time lag. For instance, if we schematize the mono- 
propellant case as a heating of the droplets from the in- 
jection temperature to a fixed deflagration tempera- 
ture followed by an instantaneous transformation in the 
final products, the time lag will be the time necessary 
to supply a definite amount of heat to the droplets; 
and if the heat transfer is taken proportional to the 
instantaneous value of the pressure to a certain ex- 
ponent n, the time lag of those particles that reach the 
deflagration temperature at a time ¢ will be given by 
the equation 


f = 7p" = [3.1] 


t’ being the integration variable. We now assume 
that the same relation applies to the whole of the proc- 
esses responsible for the time lag, so that Equation 
[3.1], or more generally 


T = To + T1; T = 70 + =7 = const; 


p™ dt’ = = const..... [3.2] 
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defines the time lag when the pressure is variable. 

If p is constant, Equations [3.1] and [3.2] are reduced 
to the first relations we have written. But this time 
even for small variations of pressure the two equations 
remain fundamentally different. In fact, for a given 
time lag, Equation [3.1 ]assumes that the variation of the 
time lag depends on the pressure values during all the 
time lag; Equation [3.2] assumes an effect of the pres- 
sure values during only one part, 71, of the total time 
lag. From the discussion of Section 2, it would seem 
that a law of the type [3.1] should represent sufficiently 
well the behavior of a monopropellant; but for the bi-- 
propellant case, the law [3.2] would be more representa- 
tive of the actual phenomenon, 7» being related to the 
mixing delay, independent of the pressure variations. 
On the other hand we have explicitly observed that 
even the pressure values during the mixing delay can 
have an effect on the rest of the time lag because they 
contribute to the combustion in “conditioning” the 
propellants. To take this effect into account a more 
complicated formula than [3.1] or [3.2] would be re- 
quired. But at this stage of the knowledge of these 
questions, it is felt that such a refinement would not 
be justified, and that only a general consideration of 
the effects of the time lag variation is required. For 
this reason, it was decided to discard even Equation 
[3.2], which involves three constants, and to select for 
the following developments the simpler Equation 
[3.1], which contains only two constants and allows 
easier mathematical treatments. 


4 Equation of the Combustion Chamber 


Before proceeding to the quantitative discussion of 
the combustion stability based on Equation [3.1], let 
us see the qualitative result of the variation of the t'me 
lag. Equation [3.1], or any other equation expressing 
the relation between time lag and pressure, shows that 
if the pressure is oscillatory, the time lag will also be 
an oscillating quantity around an average value. Now 
suppose for simplicity that the injection rate is con- 
stant. If the time lag is oscillating, it is immediately 
seen that when the time lag is decreasing the burning 
of the particles that were injected later will catch up 
with the burning of those injected earlier, and an in- 
creased rate of burning (with respect to the average) is 
the result. The opposite happens when the time lag 
is increasing. If now the increased or decreased rate 
of burning coincides with a pressure increase or de- 
crease (which happens when the average time lag and 
the period of pressure oscillations are in a certain rela- 
tion), then self-excited oscillations can exist and give 
rise to rough combustion, even in the absence of any varia- 
tion in the injection rate. 

Let us now write the equation expressing the mass 
balance in the combustion chamber. If the working 
pressure is rigorously constant and equal to /, the time 
lag will be constantly equal to 7 and the injection rate 
at every instant will equal the burning rate and the ex- 


haust rate. Let us call m this equilibrium flow rate. 
By the first approximate assumption on the time lag (Sec- 
tion 3), the propellants do not contribute to the cham- 
ber pressure until, after the time lag has elapsed, they 
are suddenly converted into gases at the final tem- 


perature. Neglecting the volume of the propellants P 
during the time lag, these gases fill all the volume of the 4 
combustion chamber. The mass M, of hot gases con- 
tained in the chamber at constant pressure p can then P 
be computed. Now suppose a deviation from the 
equilibrium conditions, and assume that the pressure 
variations are sufficiently slow so that at every instant 
the pressure is uniform throughout the combustion p 
chamber and equal to p(t). If the temperature in the [3 
chamber is supposed independent of the pressure varia- ui 
tions,‘ then the new values of the mass M, and of the - 
exhaust rate m, will smply vary proportionally with fir 
p: 
[4.1 
P ms th 
If the burning rate in the interval of time dt is called 
m,(t), the injection rate m,(t), and if the time lag is 7(t), 
the mass burning during the interval dt must be equal 
to the mass injected during the interval d(t — 7), so wh 
that of 
ty dt = d(t — with rity = = mi (t— 7)..[4.2] 
sm. 
(For brevity of writing we have here made theconvention I 
that the superscript “ means that the corresponding a 
quantity has to be evaluated at time t — 7). Now en 
the balance of mass in the chamber is 
dt = m, dt + dM, 
which can be written, using Equation [4.2]: bet) 
res 
This can be put in nondimensional form dividing by pon 
m, using Equation [4.1] and introducing the fractional E 
variations of pressure, injection rate, and burning rate ‘o 
with respect to the equilibrium values: ” 
P m m insta 
so that the mass balance equation becomes sum] 
(1 - ar) (u(t) + 1). . [4.5] 
Observe that M,/m is nothing else than the gas resi- 
dence time @, in the combustion chamber in equilibrium ane 
conditions. As the residence time before the trans- 
formation into hot gases is equal to the time lag, the Th 
total residence time of the propellants in the combus- applic 
tion chamber in steady operation is of the 
The quantity @, has an essential role in the problems of ff the ec 
combustion stability, which was disregarded by Gun- of the 
der and Friant (1), but taken into account by Yachter duce 
‘ The effects of pressure oscillations on the temperature will bes 


be considered in Section 11 to appear in the next issue. 
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(2) and especially by Summerfield (3) who used the 
expression 


6 = RT, = 


in terms of the characteristic length and velocity and 
of the temperature of the product of combustion or the 
adiabatic index y. 

6, can be used to obtain a reduced time and a reduced 
time lag: 


Coming back to Equation [4.5], let us find an ex- 
plicit expression for d7r/dt. Differentiating Equation 
[3.1] with respect to t, and assuming that p is a func- 
tion only of the integration variable t’ and not of the 
instantaneous value of 7 (an assumption to be checked 
finally), we obtain 


that is, by Equation [4.4], 


where we have again used the superscript convention 
of Equation [4.2], the last step being allowed if we 
consider only small values of ¢g, and consequently 
small relative differences between 7 and 7. 

If this expression is introduced in Equation [4.5] and 
we assume that u also is a small quantity, so that we 
can neglect the product ug, the result is the equation 


d 
g¢=m= p(d) + n(¢e = [4.9] 


between the two variables ¢ (fractional variation of 
pressure), » (fractional variation of injection rate), and 
the reduced time, z (Equation [4.7]), with the conven- 
tion ph = pw (2 — 6), and thesame for 

Equation [4.9] is the fundamental equation of the 
combustion chamber. We see now that the time lag 
appears in this equation only through 6, which is re- 
lated to the equilibrium value of the time lag; this 
means that g and therefore p are independent of the 
instantaneous values of the time lag, so that the as- 
sumption introduced in the derivation of Equation 
[4.8] is justified (this is, however, only true for small 
values of ¢). 


5 Constant Injection Rate and Intrinsic Insta- 
bility 

The simplest case to which Equation [4.9] can be 
applied is when the injection rate is made independent 
of the pressure variations. This can be obtained using 
volumetric pumps or some other device (5) provided 
the compressibility of the propellants, or the elasticity 
of the tubing, or the presence of bubbles does not intro- 
duce an appreciable variation of volume downstream 
of the device. In this case » is zero at every instant 
and Equation [4.9] is reduced to 
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contained in (6); some methods of solution are discussed in (1). 


— + =0............ [5.1] 


This is not simply an ordinary differential equation, 
since it contains the values of the variable ¢ at two 
different instants;> with the consequence that, while 
to determine the solution of a first-order linear equation 
we need only a boundary condition at z = 0 (for in- 
stance, the value of ¢); in the present case we need to 
know the values of ¢ in the complete interval of time 
fromz = —étoz=0. 

In order to find particular solutions of Equation 
[5.1], let us try the exponential ¢ = exp(az). Making 
the actual substitution we find that if the coefficient 
a satisfies the transcendental equation 


then exp(az) is a solution of Equation [5.1]. If ais a 
real quantity, the discussion of Equation [5.2] gives 
the following results. For a fixed value of n, if 6 is 
smaller than a certain 6, (which depends on n), there 
are two roots a < 0. On the other hand, if n > 1 and 
6 > than a certain 6: (depending on n), there are two 
roots a > 0. In the first case the exponential dies out 
with the time, but in the second case ¢ will increase ex- 
ponentially with time, and the solution is unstable. 
We see that the prerequisites for this type of aperiodic 
instability are n > 1, 6 > 62, so that if m < 1 (as appears 
to be likely) the exponential solutions, when they exist, 


are always stable. When n = 0, a = —1 represents the 
trivial case in which there is no influence of the time 
lag. 


Let us consider now the more interesting case of 
complex a = \ + ww. Upon substitution, Equation 
(5.2) splits into the two real equations 

A+1—n + ne- cos wd = 0) 

w — ne— sin wi = 0 \ 
Suppose we know a couple of values A, w satisfying 
Equation [5.3] for given n, 6. Then also the couple 
A, —w satisfies the equation, which means that any 
linear combination of the two conjugate quantities 
exp(az) and exp(&z), or—remaining in the real field— 
any linear combination of e“ sin (wz) and e” cos (wz) 
satisfies the equation. This combination represents a 
real oscillatory solution ¢ with arbitrary phase, oscil- 
lating with a reduced angular frequency w with respect 
to the reduced time z, or an angular frequency Q = 
w/6, with respect tot. The period is 7 = 2r/Q2 and the 
frequency, the inverse of this quantity. The quantity 
wd which appears in Equation [5.3] has the following 
simple interpretation; because of Equation [4.7]: 


so that wé, measured in revolutions, represents the ratio 
between the time lag and the period of the oscillations. 


5 This kind of equation is sometimes designated as hystero- 
differential equation. The discussion of a particular case is 
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@ The amplitudé of the oscillation varies with time fol- 
lowing the factor exp(Az); if \ < 0, the amplitude de- 
creases and the solution is stable; if \ > 0, the ampli- 
tude increases exponentially and the solution is un- 
stable. The neutral case, \ = 0 represents, therefore, 
the limiting case between stable and unstable solutions, 
when the amplitude remains constant. 
From the two Equations [5.3], we find 


wd 
64+ (1 —n)ji= “fan 
Introducing the value of \é obtained from this into the 
second Equation [5.3], we have the following equa- 
tion: 


wd 
6 — 
= = tan — 
f(6, n) = née = h(wd) 


which can be solved graphically, for any n, 6, plotting 
h(wé), as we have done on the right side of Fig. 1, and 


the first, or fundamental, mode; but more than one 
for all the other higher frequency modes. For small 
5 or large n, the first root can become imaginary, but 
the higher frequency modes exist always. 

Considering first the case n = 1, we see that as f(1, 6) 
coincides with né = 6 the solution will be stable, neu- 
tral or unstable if k 2 h = 6. The solution is therefore 
stable when wé is conta‘ned in the first quadrant, un- 
stable when it is in the second quadrant. Hence, the 


fundamental mode is stable or unstable if 6S 6% = Z 
the first higher frequency mode is stable or unstable if 
6$ 5 a etc. When 6 increases from small values the 


first solution to become unstable is the fundamental 
mode. When this one is unstable, the fact that the 
others are still stable or already unstable has but little 
importance. We conclude that the examination of the 


| 


kwd) 


| 


\ hiwd) 


Fig. 1. 


cutting the curves with a horizontal line of ordinate 
f(6, n). This quantity is shown as a function of 6 on 
the left side of Fig. 1 for three different values of n. The 
value of \ can then be found using, for instance, the 
second Equation [5.3]: 


sin wd nb 


AS = 


The function k(w6) is plotted on the right-hand side of 
Fig. 1, and né is shown on the left-hand side. We see 
that the solution is stable, neutral, or unstable if nd $ 
k(wd). The examination of Fig. 1 shows that for every 
couple of values n, 6 there is an infinite number of roots 
w6, corresponding to an infinite series of increasing w, 
the first root being contained between 0 and 7z, the 
second between 2z and 3z, etc. There is less than one- 
half oscillation in the duration of the time lag for 
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GRAPHICAL SOLUTION OF EQUATION [5.3] 


behavior of the solutions for 0 < wé < 7 gives in this 
case all the required information about stability. 

The case n < 1 is illustrated in Fig. 1 for n = 0.75. 
In this case, né < f(n, 6), so that the stability can exist 
also when h > k, and the stability conditions are im- 
proved. Here again, for 6 small, all the modes are 
stable; after a certain critical 5% larger than for n = 1, 
the fundamental mode becomes unstable; and _ with 
increasing 6 other successive modes become unstable: 


and the behavior for n not too small, is analogous to 
that previously discussed. The quantitative evalua- 
tion of 5» will be deduced in the following. 

Finally, when n > 1 (as in the case n = 2 of Fig. 1) 
né > f(n, 6) and instability can exist with h < k; which 
means that the range of stability is reduced, and 6» de- 
creases with increasing n. 
clusions are the same as before. 
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A general conclusion from this discussion is that 
there are always an infinite number of solutions of 
Equation [5.3] and therefore of Equation [5.1]. A 
linear combination of these solutions, with the even- 
tual use of the exponential solutions corresponding to 
Equation [5.2] with real a, solves the problem of find- 
ing the solution of Equation [5.1] with a quite general 
shape of the prescribed initial distribution of g(z) be- 
tween —6 and 0; and the corresponding solution is 
stable or unstable following if the reduced time lag 6 is 
smaller or larger than a certain critical value 6x. 

The kind of instability discussed in this section can 
be denominated “‘intrinsic” since it depends only on 
the conditions inside the rocket combustion chamber, 
and not on variations of the injection rate. The value 
of 6x and of the corresponding ws can be deduced di- 


rectly from Equation [5.3], putting \ = 0. We ob- 
tain: 
COS (wy 5,4) = = n sin (wy | 
= Vn? — (1 —n)? = V2n—-1; | 


= 


The corresponding values are plotted against n in Fig. 
2 

Real values are obtained only when n > 0.5, which 
means that below 0.5 the solutions are always stable. 
The intrinsic instability can therefore exist only if n > 
0.5. 
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Let us finally observe that if, instead of Equation 
[3.1] we select Equation [3.2] to represent the pressure 
dependence of the time lag, the results are the same 
provided m replaces n and 7; replaces 7, that part 79 of 
the total time lag which does not vary with pressure 
being in this case without influence on the stability. 


6 The Feeding System 


The assumption of constant injection rate used in 
the preceding section does not correspond generally to 
the real operating conditions of a rocket motor. The 
oscillations of the chamber pressure are in fact respon- 
sible, with the ordinary feeding systems, for oscilla- 
tions in the injection rate which have in turn an effect 
on the pressure oscillations. This may be true even 
when using devices designed to produce a constant mass 
flow, placed at a certain section of the tubing; since 
variations of the mass flow may still occur between 
this section and the injector if some intermediary part 
has an elastic behavior (gas or vapor bubbles, compres- 
sibility of the propel'ant, elasticity of the walls). We 
exclude from our analysis the case of a feeding system 
with centrifugal pumps, whose response to chamber 
pressure variations cannot be defined in a general way, 
especially when some kind of automatic control is used. 
In the last hypothesis, oscillations in the automatic 
control system can arise that may interfere with the 
operation of the combustion chamber, in which case 
the analysis becomes very involved. 

To confine our considerations to simple and general 
cases, let us consider two types of feeding system: (a) 
A constant rate supply system, connected with the in- 
jector by lines of a certain length and with an elastic 
behavior; and (b) a constant pressure supply system, 
again connected with the injector by elastic lines. For 
the moment, we will consider a monopropellant system 
to avoid complications. In Section 10 (Part 2) we 
shall consider separately the effects of a bipropellant 
combination. 

Case (a): Fig. 3 shows this type of feeding system 
schematically. When the flow is steady, neglecting the 


friction and other eventual losses,® we have 


CONNECTING LINE| 


SCHEMATIC REPRESENTATION OF A CONSTANT RATE FEED- 
ING SYSTEM WITH ELASTICITY IN THE LINES 


FIG. 3. 


6 These can be taken implicitly into account in this case and 
in some other cases by using instead of A; an equivalent area, 
such that the Equation [6.1] is still satisfied. 
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where A, represents the area of the jets from the in- 
jector, that is, the orifice area corrected for the con- 
traction coefficient. Under the pressure p; the capaci- 
tance, C, takes a constant volume so that the mass flow is 
the same before and after the capacitance. When the 
flow is unsteady the pressure will vary and also the 
capacitance C; for small variations of pressure C will 
vary linearly with p,, and we can define a coefficient 


In the supply line, the mass flow is constantly m; the 
continuity equation then gives for the mass flow in the 
connecting line (injection rate) : 


Equating the work performed per unit time by the in- 
stantaneous pressure drop p, — p to the rate of in- 
crease of the kinetic energy in the connecting line, plus 
the kinetic energy produced in the jets per unit time 
(neglecting the small variations of kinetic energy be- 
tween the supply line and the connecting line, and in 
the capacitance itself) we obtain: 


(pi — p) = + [6.4] 
with 

nm = pAw = pA; M; = plA:; ce [6.5] 


where M, represents the fluid mass in the connecting 
line.’ 
Equation [6.4] can be written 


[6.6] 


Subtracting Equation [6.1] from Equation [6.6] we ob- 
tain: 


Equation [4.4], and the three quantities: 


the preceding equation becomes: 


7 When the connecting line is composed of pieces of different 
length and areas, and we select a representative area A;, with the 


corresponding velocity v:, M; is defined by: M,.(v;2/2) = es Mi 
(v,2/2); andl by the second Equation [6.5]. 


the last step being justified for small variations of mass 
flow. Dividing by 2Ap; introducing the reduced time 
z (Equation [4.7]), the fractional variations defined by 


Using the same variables, Equation [6.3] becomes 
dy 
p= -E [6.9] 
where the elasticity parameter 


[6.10] 
mh 8, 
includes all the effects of the elasticity of the system, 
and can be interpreted as the double of the ratio be- 
tween the increase in volume of the system under the 
pressure Ap and the volume injected during the gas 
residence time. 

Equations [6.8] and [6.9] are sufficient to describe 
the behavior of the feeding system in case (a). They 
depend only on the three constants, P, J, E. The 
meaning of E has just been discussed; the pressure 
drop parameter P has an immediate interpretation as a 
measure of the relative value of the pressure drop 
through the injector; and the inertia parameter J, 
which includes all the effects of the inertia of the pro- 
pellant in the lines, can also be written: 


so that it represents the ratio between the kinetic energy 
of the propellant in the connecting line and the work 
performed by the pressure drop during the gas residence 
time, or (what is the same) the kinetic energy produced 
in the injected propellant during the gas residence time; 
everything being computed under steady conditions. 
This definition of J remains true for composite connect- 
ing lines. 

Case (b): This type is shown schematically in Fig. 
4. Again neglecting the friction in the lines and the 
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FIG. 4. 


FEEDING SYSTEM WITH ELASTICITY IN THE LINES 


and in nonsteady flow the two equations, corresponding 


to Equation [6.6]: 


SCHEMATIC REPRESENTATION OF A CONSTANT PRESSURE 


variation of kinetic energy from the supply to the con- 
necting line, as well as the kinetic energies in the ca- 
pacitance and in the tank, we have in steady conditions: 


m2 
Po — Pi = 0; Po — P = = Ap. 
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P A, dt’ 


. [6.12] 


where m, and m represent the mass flows in the supply 
and in the connecting line, whose difference satisfies to: 


d 


similar to Equation [6.3]; and where the pressure po 
generated by a regulated gas pressure is assumed con- 
stant. 

Subtracting Equations [6.11] from [6.12], dividing by 
2Ap, introducing the quantities defined by Equations 
[6.7] and [6.10] and the fractional variation of mass 
flow from the tank 


Equations [6.12] and [6.13] are transformed into 


dys, 
y= “2; 
d 
Eo 


which are sufficient to describe the behavior of the feed- 
ing system in case (b), and again contain only the three 
constants P, J, E whose meaning is the same as in 
case (a); plus the supplementary parameter £. 


7 Limits of Stability with Constant Rate Supply 
(Case (a)) 


iquations [4.9], [6.8], and [6.9] are sufficient for the 
description of the behavior of the three unknown func- 
tions g(z), u(z), ¥(z) when proper initial conditions are 
given. In fact, they could be reduced by elimination to 
a single third-order hystero-differential equation in, say, 
¥(z) alone, an equation that could be treated in a way 
similar to that of Section 5. Instead of actually carry- 
ing out the elimination, we can write the equation 
corresponding to Equation [5.2] by trying directly in 
the system a particular solution of the type: 


a, A, B, C being complex constants. Eliminating the 
common factor exp(az), we see that Equation [7.1] 
represents a solution of our system if the three homo- 
geneous, first-degree equations in A, B, C 
(a +1 —n + ne-25)A — =0 
PA+(1+Ja)B-—C=0 
B+ EaC =0 
are simultaneously satisfied, which happens if, and only 
if, the determinant of the coefficients vanishes, that is, 
if 


—e —ai 0 
P 1+Ja —-1/=0 
0 1 a 


Expanding the determinant we find that if « satisfies 
the equation 
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(1 + Ea + JEa?) (a +1 — n + ne- + 9) 
then Equation [7.1] represents a possible solution of the 
system. Equation [7.2] is reduced to Equation [5.2] 
when E = 0, and in general it can be discussed in the 
same way; also in this case real roots can exist, but the 
most interesting are the complex conjugate roots a = 
\ + iw, which correspond to oscillatory solutions, stable 
or unstable following if \ $ 0. Again, an infinite num- 
ber of roots are possible, corresponding to increasing 
frequencies and dampings; and again the solution cor- 
responding to the first root will be the first to become 
unstable when 6 increases from small values and reaches 
a value és, w taking a certain value wx. The deter- 
mination of these critical values can be performed di- 
rectly from Equation [7.2], upon substitution of a = 
We obtain: 
(l—n-+ (1 — JEw,? + = 

—[n(1 — JEws?) + i(n + P) inde 
Equating the moduli and the arguments of the two sides. 
of this equation and considering that for the first root 
w, 6, must be <27, we have 
[((1 — n)? + we?) — JE + = 

— JEws?)? + (n + 


E. 
x + tan“ 


The first of these equations is a third-degree equation 
in wx?; if wx? is deduced from here in terms of n, P, J, E, 
and introduced in the second Equation [7.3], this would 
give the value of 6x in terms of the same four quantities, 
and we would obtain the critical value of the time lag 
and of the angular frequency for any given combi- 
nation of the constants. But the solution of the system 
is simpler in the following way. Solve the first Equation 
[7.3] for 1 — JEws?: 


(n + P)? — — (1 — _ 
+ 


1 — JEw,? = 
. [7.4] 
(® representing for brevity the inverse of the radical), 


and replace this value in the second Equation [7.3]; 
solving for 6x we obtain: 


= —| — 


which shows that the relation between 6x and wx de- 
pends only on the values of n and P, but not on J and 
E. Figs. 5, 6, 7 show this relation for n = 0, 0.5, and 
18 computed from Equation [7.5] for various values of 
Ap/p = 1/2P. On each curve there are two values of 
the ordinate for a given value of w*, corresponding to 
the two signs of Equation [7.5]. The reality of the 
radical @ involves the inequality 

8 These values have been selected only as examples; the ac- 


tual values of n have to be derived experimentally, and might 
be even larger than unity. 
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2n — 1 = n? — (1 — n)? < we? < (mn + P)? — (1 — nn)? 

= (P + 1)(2n + P - 1) 
so that the curves are all contained between a mini- 
mum and a maximum of the abscissa. The minimum 
is different from zero only when n > 0.5; when n < 0.5 
the curves extend up to the ordinate axis where the 
critical time lag goes to infinity. In the last case, a 
limitation exists also for P since, in order to have a 


positive value of the maximum w,’, the inequality P > 
1 — 2n must be satisfied. When, on the contrary, 
P < 1 — 2n, no real value of ws and 6« can be found, 
and the stability of the combustion is assured no mat- 
ter what values 6, E, J take. When n = 0 this condi- 
tion coincides with the result of (1) and (3) which show 
that for their particular feeding system, the com- 
bustion is always stable when Ap/p > 0.5. We see now 
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that the possibility of assuring an unconditioned sta- 
bility through an increase in Ap exists only when n < 
0.5; and that the necessary increase in Ap becomes more 
and more important when n approaches 0.5; until for 
n > 0.5 even an infinite increase of Ap may not be 
sufficient to produce stability. Remembering the re- 
sults of Section 5, we recognize that this happens when 
the operation can be intrinsically unstable, in which 
case it is not possible to obtain stability by changes 
in the parameters of the feeding system. 

In Figs. 8, 9, and 10 we have plotted for n = 0, 0.5, 
1 the values of wxdx/2r = t+*/Tx against 5x. Here 
again the curves Ap/p = constant are closed curves if 
n > 0.5; while they extend to infinite values of 6. if 
n < 0.5. When n > 0.5 the points corresponding to 
the minimum value of ws (which makes = 0) have 
the co-ordinates 


— 
n—l1 


V2n 
independent of P, and all the curves pass through the 
same point. When n = 0.5 this point goes to 6x = ©, 
and becomes imaginary when n < 0.5. 

For all values of n the co-ordinates of the points cor- 
responding to the maximum value of w*« (which makes 
o)are 


73 l-—n 1; 
tant 4/(" +?) —1 


V(n +P)? — (1 — 


(54 


(64) 
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The curve connecting these points for various P sepa- 
rates the values corresponding to the upper sign of 
Equations [7.4], [7.5] (above the curve) from those cor- 
responding to the lower sign (below the curve). For 
n = 1 the curve becomes an horizontal line. 

At every point of a given curve Ap/p = const, 
wx and ©® take definite values. The selected point 
represents the critical conditions with a given feeding 
system provided the values of E and J of this system 
satisfy Equation [7.4]. The simplest case is when the 
length of the connecting line is negligible and J = 0. 
Then, from Equation [7.4], 


Wx 


But by its definition (Equation [6.10]), Z must be posi- 
tive, so that in this case only the upper portions of the 
curves have a physical meaning. The point ® = 0 
corresponds therefore to Ey) = 0, and the line 6 = @ 


to Fy = «. Lines of equal Zo have been drawn in 
Figs. 8, 9, and 10. When n < 0.5 the point @ = 0 is 
imaginary, which means that at 2) = 0 the conditions 
are always stable (in accordance with the results of 
Section 5). In the case n < 0.5, Equation [7.6] shows 
that Eo = © both when @ = © and when ws = 0 
(6% = ©) so that if we move in the clockwise direction 
on the curves of, for instance, Fig. 8 starting from the 
line Hy = ©, Ep will first decrease, reach a minimum 
value Eomin, and then increase again. The combus- 
tion is stable for all 6’s when Ey < Eomin} 9S SOON as 
Ey > Eomin the combustion becomes suddenly unstable 
if 6 is larger than the value of 6x corresponding to Eon} 
for Ey > Eomin there are on the curves two points cor- 
responding to the same value of Ey. The one contained 
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in the portion of the curve between Ey) = © and Ey = 
Eomin has the smaller 6+ and therefore is the first to be- 
come unstable, and the only one of practical importance. 
When n increases the value of 5% at which Eo = Eoin 
increases, until it reaches infinity when n = 0.5, still 


with a finite value of #) = 1/V P(P + 1); as soon as 
n > 0.5 the curves become closed and at the value of 
6x corresponding to ® = 0 (which decreases from in- 
finity with increasing n) we have Ey) = 0. For n > 
0.5, therefore, instability exists even with Ey = 0 if 
6 is larger than the corresponding 6x. From Fig. 10 
we see that in this case for small values of Eo, 6« in- 
creases slightly and the stability is improved; a fur- 
ther increase of Eo, however, makes 6% drop again, and 
the system becomes less stable. This is true for all 
the values n > 0.5. 

In the general case J # 0, Equation [7.4] can be 
written: 

= + J Wy 

Ey) now being only the mathematical expression de- 
fined by Equation [7.6], without any physical corre- 
spondence, so that also negative values are allowed. 
The point E = 0 coincides with Ey) = 0 (@ = 0), when 
this point exists, that is when n > 0.5; on the other 
hand, when n < 0.5 and wx = O(é% = ©), E is still 
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infinity like Zo, no matter what is the value of J. How- 
ever, E = = is also obtained, for all values of n, when 


and as J must be positive, it is obtained for negative 
Eo, below the line Ey = ~. For n < 0.5, between the 
two infinite values of E there is a minimum, smaller 
than E ,,;, but still positive, and again FE = 0 gives 
certain stability. When £ is larger than the minimum 
there are two corresponding values of 6%, of which only 
the smaller has an actual interest. In the correspond- 
ing region when n < 0.5, and everywhere when n > 0.5, 
the effect of an increase of J is to displace the critical 
point in the counterclockwise direction on the curves 
of Figs. 8, 9, and 10; and therefore an increase or a de- 
crease of 6s (with a better or worse condition of sta- 
bility) may follow, depending on the location of the 
point on the curve. However, the maximum and the 
minimum possible values of 6 remain unchanged, and 
depend only onnand P. ~ 


8 Limits of Stability with Constant Pressure 
Supply (Case (b)) 


This case is described by the four Equations [4.9], 
[6.14]. Here again, a process of elimination would 
lead to a single fourth-order hystero-differential equa. 
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tion in, say, mi(z), which is reduced to a third-order 
equation when € = 1 or when — = 0. If the concen- 
trated variable capacitance of Fig. 4 is assumed to repre- 
sent the diffused elasticity of the line and the enclosed 
propellant, the proper location of it would be in the 
center of the line, that is, § = 0.5. However, we have 
not developed this case, and we have discussed only the 
case £ = 1 because the corresponding formulas can be 
simply discussed using the treatment of Section 7 (the 
case £ = 0 is a trivial case in which the elasticity has no 
effect at all on the system; and can be included in the 
f= 1 case assuming HF = 0). 

Using, with € = 1, the same technique of Section 7, 
we find that if the four unknown functions are assumed 
to vary proportionally to exp(az), a must be such that 
the determinant 


ja +1 —n + ne-ai —e—ad 0 0 | 
1 0 —1 
0 1 -1 Ea 
0 0 Ja 1 


has zero value. 
Developing the determinant we obtain the equation 


(l + Ja + JEo?) (a +1—n + ne~%) + 
PUL + = 0 


whose discussion is conducted in the same fashion as 
for Equations [5.2] and [7.2], and gives analogous re- 
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sults. If we are interested only in the limit of stability 
(when the fundamental mode becomes unstable) we 
substitute again a = iw*, and the preceding equation 
gives 
(1 —n+ iw,) (1 — JEo,? + = 

—[(n + P) (1 — JEa,?) + inJw,y]e— [8.1] 
Equating the moduli and the arguments we obtain: 


+ [(1 — JE + = 
(n + P)*1 — + 


.. [8.2] 
J ws 


n 
n+P1 — JEw,? 
From the first Equation [8.2] we have 


a + 


+ (1 — — 

(n+ P)?—w,? — (1 — n)? 

Ja,®. eee [8.3] 

® being the same function as in Section 7; replacing in 
the second Equation [8.2]; and solving for dx: 


1 n 


Wx % 


1 — JEw,? = F Joy 


which is exactly the same as Equation [7.5]. We con- 
clude that as far as the relation between 64 and w« for 
given values of n, P is concerned, the system of Fig. 4 
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with = 1 behaves exactly like the system of Fig. 3. 
The curves of Figs. 5, 6, 7, 8, 9, and 10 represent, there- 
fore, also the present case, abstraction being made 
from the actual values of Z and J. A given point on 
one of the curves represents the critical condition for a 
certain system if the relative E, J satisfy to Equation 


[8.3]. The simplest case when the elasticity is neg- 
ligible, E = 0, gives 


As in this case Jo represents the essentially positive 
physical quantity defined by the second Equation [6.12], 
only the lower sign has a physical meaning, and there- 
fore only the lower portion of the curves (below the 
@ = © line) corresponds to a real process. Lines of 
equal Jp have been drawn in Figs. 8,9, and 10. At the 
point = 0 we have Jp = ~; and the line 6 = © 
corresponds to Jo = 0 which coincides therefore with 
the Ey) = © line. For given n, P, small values of 
Jo (moderate length of lines) result in a decrease of 5 
and of stability; a further increase of Jo up to Jo = © 
results in an increase of 6s and an improvement of sta- 
bility. When n < 0.5 above a certain J» the com- 
bustion is always stable; but if n > 0.5 even when Jp = 
©, 6x is finite and the combustion can be unstable. 
An increase of Ap/P at constant Jo results generally in an 
increase of 6« and an improvement of stability, but the 
combustion will remain unstable even with Ap = & 
when the rocket is intrinsically unstable. 

The general case E ¥ 0 can be easily discussed ob- 
serving that Equation [8.3] can be written 


= + Ew,? 
where positive or negative values of Jo are possible, 
since J) now loses its physical meaning and is only the 
mathematical expression of Equation [8.4]; we see 
that at Jo = 0, J is always zero; and J = © when 

E, being the function defined by Equation [7.6]. Hence 
the lines of constant Ey of Case (a) correspond also 
to the conditions E = Ey, J = @ of Case (b). (Simi- 
larly, the J defined by Equation [7.7] coincides with 
J, defined by Equation [8.4]; so that the Jo = constant 
lines of Case (b) represent also the conditions J = Jo, 
E = o of Case (a).) It follows that the result of an 
increase of E is a counterclockwise displacement of the 
critical point on the corresponding curve of the figures, 
with a resulting increase or decrease of 6x, depending 
on the location of the critical point. If this point is 
in the most probable region of J (that is, between 6+ min 
and 6% max) the added elasticity will produce an increase 
in 6x with beneficial effects on the stability. For given 
n, P and for every J there exists an optimum value of 
the elasticity which renders 6s maximum. 
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9 Effect of the Nonuniformity of the Time Lag 


In Section 3 we made the explicit assumption that 
the time lag is the same for all the particles. Let us 
now assume more generally that different fractions of 
propellant injected at the same time have different 
time lags, so that this time lag is not only a function of 
the pressure, but also of the location of the considered 
fraction in the jet. The mathematical formulation of 
this assumption can be derived simply in the following 
way. Suppose that we inject the propellant through 
N distinct injectors, each one having a uniform time 
lag, different from the time lag of the other injectors. 
If 7, is the time lag of the Ath injector (h = 1,2,...,N), 
and we assume a relation analogous to Equation [3.1] 
with the same exponent n for all the injectors, we ce- 
duce the equation similar to Equation [4.8]: 


dr, 


Now the burning rate from this injector, at time ¢, is 
related to the injection rate at the time t — 7» by the 
equation similar to Equation [4.2]: 


d 


and the total burning rate at time ¢ will be 


rin PLL ne — 


If we pass to the limit of an infinite number of injectors 
and we replace min,“ with m°“ df where: f, variable 
between 0 and 1, represents the location of the fraction 
we are considering in the jet, 7 is now a function of f, 
and m;‘") represents the total injection rate at the 
time t — r(f); and if instead of the summation we per- 
form an integration, we can write: 


1 


Let us now introduce the fractional variations of the 
burning and of the injection rate defined by Equations 
[4.4]; we have, neglecting the terms of higher order 
(products of and ¢): 


1 
ng +f, || | [9.1] 


where we have already replaced the time lag with the 
(nonuniform) reduced time lag 6(f) (Equation [4.7]). 
Introducing Equation [9.1] in the first two members of 
Equation [4.5], we arrive at the equation of the com- 
bustion chamber: 


dz 0 


which for uniform 6(f) becomes Equation [4.9]. 

If we try on this equation a solution of the type of 
Equation [7.1] we find, cancelling the common factor 
exp(az): 
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Equation [9.2] replaces, in the case of nonuniform time 
lag, the analogous equation derived from Equation 
[4.9] and used in Sections 7 and 8. We see that from 
the equations developed there we can obtain those rela- 
tive to the nonuniform time lag case by simply replac- 
ing everywhere exp(— a6) with the integral 


which depends, of course, not only on a but also on the 
distribution function 6(f) of the time lag. 

To discuss in some simple cases the effects of this 
substitution, let us define an average value of the re- 
duced time lag 


1 
in = 
and the variation from the average 
1 
df) = WS) Af [9.4] 
The integral [9.3] then becomes 
Ia) = fy df = Be-aim.... [9.5] 


where the coefficient 8(a) defined in this way will now 
contain all the effects of the nonuniformity of the time 
lag when we compare cases with the same 6,,, and is 
unity for uniform 6 = 6». 

If we are interested in finding the limit of stability 
of a given system we can replace a = iwx and obtain 


B= cos (wye) df — if, sin (wye) df....... [9.6] 


Suppose now that f is arranged in such a way that 
i(f), and therefore ¢(f), increase monotonically with 
f. Then, due to Equation [9.4], from a certain nega- 
tive value « at f =0, eincreases to a positive value eat 
f = 1. If this curve is distributed in an antisymmet- 
rical way with respect to f = 0.5, that is, if e(f) = 
—e(1 — f), then the second integral of Equation 
[9.6] is zero, and @ is real and smaller than one. If 
the curve is not antisymmetrical there will be an imagi- 
nary part of 8, but probably for plausible ¢(f) it will 
be small with respect to the real part, still <1. The 
actual computation of 6 requires the knowledge of 
«(f). Of course this function depends on the particular 
injection system and on the type or combination of 
propellants; but it is probable that the variation of 
§ and ¢ is not uniformly distributed with f (linear vari- 
ation), and that a larger fraction of propellant has 6- 
values close to 6m, or e-values close to zero, and a smaller 
fraction has e-values close to the two extremes €, €1. 
In order to evaluate numerical examples, let us assume 
two different e(f), both with an antisymmetrical varia- 
tion so that = — «1: 
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the second assumption being a simple example of the 
distributions just discussed. The value of 8 can be 
readily evaluated in both cases and is 


cos Wy €1 
2 
(2 
where ¢; = (6; — 69)/2 is half the total excursion of the 
reduced time lag. It is again checked that B is < 1. 


The minimum value is obtained when e« is maximum, 
that is, when 6) = O and e = 6;/2 = 6», so that 


Wy €1 


COS bm 
Bmin = 2 
1 


We can use these formulas in the simplest case of the 
intrinsic instability (constant rate of injection). The 
corresponding equations are obtained from Equation 
[5.3] by introducing 6,, instead of 6, and 6 exp(—A6é,,) 
instead of exp(—)é) and finally putting \ = 0. We 
obtain 


sin Il: 


Wx bm 


... [9.8] 


l-—n 
B COS wybmy = — 


‘ 
B sin = 


If n = 1 the first always gives w*6,* = 2/2, and the 
second ws = B, so that 


bmx 


2B 


nla 


and the critical time lag (average) is increased with 
respect to the case of uniform 6. The amount of the 
increase depends on the value of 8. If we take the 


minimum possible values of Equations [9.8] (4 = 
5m), We have 


us 


I: Bmin= — = 0.637; II: Bmin = an 0.786 


us 
and the corresponding increases of 5*« are of 57 per 
cent and of 27 per cent. 

The value of « = 6, just used is probably larger 
than the values likely to appear in practice. As a 
rough estimate, we can say that if the total residence 
time of the propellants is computed from Equation 
[4.6] with the average value of 7, the total reduced resi- 
dence time is 


im + 1; 

on the other hand, with « = 6,,, the maximum reduced 
time lag is 26,, which, with 6,, = 1/28, is larger than 
the total reduced residence time, and the correspond- 
ing fractions of propellant would not have the time nec- 
essary to burn; that is, a loss in combustion effi- 
ciency would result. If we assume a more conserva- 
tive value, «, = 6,,/2, which does not produce the 
effect just discussed, we have from Equation [9.7]: 


2 3 


and the increases of 6,,*« are of 11 per cent and 6 per 
cent. 
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We conclude that the effect of the nonuniformity of 
the time lag is in this case an improvement of the sta- 
bility, and that the improvement for distributions of 
the more probable type II and not too extreme excur- 
sions of the time lag is only of a few per cent. 

Analogous conclusions can be reached for other 
values of n. For instance, the minimum value of n 
for the existence of intrinsic instability is raised with 
respect to the uniform case. 

From these simple computations, and from others 
that we exclude for brevity, it seems possible to con- 
clude that the nonuniformity of time lag has always a 
favorable, though not too important, effect on sta- 
bility. 
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cooled rocket motors. 


Address inquiries to: 


Special Refractories 
CARBORUNDUM 


Special silicon carbide refractories made by The 
Carborundum Company have proved to be the 
best available materials for ceramic linings in un- 


The Carborundum Company 
Refractories Division 
Perth Amboy, New Jersey 
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Chemical Kinetics and Rocket Nozzle Design 


By F. J. KRIEGER! 


The Rand Corporation, Santa Monica, Calif. 


The effect of chemical kinetics on rocket nozzle design 
was investigated for the particular case of hydrogen gas 
flowing adiabatically through a typical rocket nozzle hav- 
ing a chamber-to-throat-area ratio of 2 to 1. For com- 
parative purposes the following types of flow were con- 
sidered: (a) Constant composition (frozen equilibrium), 
(b) instantaneous chemical equilibrium (shifting equilib- 
rium), and (c) kinetic chemical equilibrium. A _ step- 
wise iteration process was employed to perform the inte- 
gration of the differential equation in case (c). In each 
case the gas entered the nozzle at a temperature of 3500 K 
and a pressure of 20 atm and was allowed to expand isen- 
tropically to an exhaust pressure of 1 atm. The results 
for a mass flow rate of 1000 gm sec”! are presented in Figs. 
2through6. The instantaneous equilibrium flow assump- 
tion not only gives a higher specific impulse, but also re- 
quires a larger nozzle than either the kinetic or the con- 
stant composition flow assumption. The kinetic equilib- 
rium flow results are intermediate between those for in- 
stantaneous equilibrium flow and those for constant com- 
position flow, the relative position depending on the mag- 
nitude of the reaction rate which governs the kinetic 
equilibrium. 


T HAS been standard practice in the field of rocketry 
to compute the specific impulse (i.e., the jet velocity 
divided by the acceleration of gravity) of a propellant 
system on the basis of either of two fundamental as- 
sumptions, namely, that the adiabatic flow of propellant 
gases from the combustion chamber through the nozzle 
is such that either (a) constant composition or (b) in- 
stantaneous chemical equilibrium is maintained. For 
an adiabatic process, assumption (a), which says that 
the inherent reaction rate is zero, gives a minimum 
value for the specific impulse, whereas assumption (b), 
which says that the inherent reaction rate is infinite, 
gives a maximum value for the specific impulse. The 
introduction of reaction rate equations into the com- 
putation of specific impulse so complicates the problem 
that heretofore investigators have, perforce, been gen- 
erally satisfied with a knowledge of the limiting values 
obtained on the basis of the above assumptions. 

Aside from the computational difficulties which it 
entails, the kinetic approach, because it requires close 
adherence to physical reality, tends not only to estab- 
lish the true specific impulse for a propellant system, 
but also to give some insight into the design of the 
nozzle for that system. In general, specific impulse 
values computed under assumptions (a) and (b) are 
made without regard to nozzle configuration and with 
the added simplifying assumption that the chamber 
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gases enter the nozzle with zero velocity. Since this 
latter assumption has no physical justification and is 
erroneous, it leads to spurious values of the specific im- 
pulse. In order to obviate this condition it is neces- 
sary at the outset to determine the chamber velocity at 
the entrance to the nozzle for the type of gas flow under 
consideration. This may be done only by considera- 
tion of nozzle configuration. 

The nozzle generally accepted by rocket designers is 
the convergent-divergent De Laval type, which is 
formed from two intersecting right circular cones. 
The half angle of the convergent cone is usually 30° and 
that of the divergent cone, 15°. To prevent discon- 
tinuity of flow, the throat section is rounded off with a 
radius of curvature which is equal to the throat diame- 
ter 2r,. Such a nozzle was considered by Altman and 
Penner (1), (2).2. The nozzle considered in this study 
is shown in Fig. 1. It differs from the Altman and Pen- 
ner nozzle in two respects. First, according to current 
design practice, the entrance area A, is twice the 
throat area A, This restriction permits one to es- 
tablish, among other things, the origin of a rectangular 
coordinate system at the entrance to the nozzle. Sec- 
ondly, the convergent cone is replaced with a sphere 
of radius r,. The sphere is tangent to the rounded 
throat section at a point where the slope makes an angle 
of 28° 31’ with the nozzle axis. The use of a spherical 


section instead of a convergent cone insures conti- 
nuity of flow from the combustion chamber into the 
nozzle entrance. 


FIG. 1 CROSS SECTION THROUGH ROCKET NOZZLE 


2 Numbers in parentheses refer to Bibliography at end of 
paper. 
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The propellant system in this study is pure hydrogen 
having a chamber temperature of 3500 K and a chamber 
pressure of 20 atm. Its history prior to its arrival at 
the entrance to the nozzle is irrelevant. From the 
basic nozzle design parameters and mass flow rate con- 
siderations, the following properties may be computed: 
(1) Chamber velocity v,; (2) nozzle throat area A,; 
(3) nozzle length L; and (4) jet velocity v, or equiva- 
lent specific impulse J. In this study, the foregoing 
properties are determined not only for (a) constant 
composition flow and (b) instantaneous chemical equi- 
librium flow, but also for (c) kinetic chemical 
equilibrium flow. Case (c) makes use of the expres- 
sion for the rate of recombination of hydrogen atoms 
during the adiabatic expansion process of the propellant 
gas as it flows through the nozzle. 


General Thermodynamic Equations 


Several basic equations may be derived from a ther- 
modynamic consideration of the propellant gas mix- 
ture. For the association reaction 


2H = 


the equilibrium constant K, is given by the expression 


where Ny, and Ny are the mole fractions of H2 and 
H, respectively, and P is the total pressure in atmos- 
pheres. The mole fractions satisfy the obvious rela- 
tion 


The molecular weight of the gas mixture is given by 
the relationship 


M => NiM: 


where N; and M; are the mole fraction and the molecu- 
lar weight, respectively, of component 7. 

The expansion of the propellant gas is assumed to be 
isentropic. The specific entropy of the gas mixture is 
given by the expression 


= (DNS! — RY Rin P)......141 


where s is the entropy in calories per degree per gram, 
M is the molecular weight of the gas mixture, S;° is the 
entropy per mole of component 7 at one atmosphere 
and at a given temperature, N; is the mole fraction 
of component 7, R is the gas constant, and P is the given 
pressure in atmospheres. 

The specific enthalpy of the gas mixture at any point 
in the expansion process is given by the expression 


where h is the enthalpy in calories per gram, M is the 
molecular weight of the gas mixture, and N;, AHy,, 
and AH7..16, ; are the mole fraction, the heat of for- 
mation, and the sensible heat content, respectively, of 
component 7. 
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By equating the change in kinetic energy of the ex- 
panding gas to its change in enthalpy, the velocity at 
any point in the expansion is found to be 


(6) 


where v is the velocity in centimeters per second, h is the 
specificenthalpy in calories per gram, J is the mechanical 
equivalent of heat, and the subscript c refers to the 
chamber conditions. 

Specific impulse and velocity are related by means 
of the simple expression 


where J is the specific impulse in seconds, v is the ve- 
locity in centimeters per second, and g is the accelera- 
tion of gravity. This concept of specific impulse as- 
sumes complete or perfect expansion of the propellant 
gas. 


Mass Flow Equations 


A fundamental relationship between nozzle con- 
figuration and propellant gas mass flow may be derived 
from the following considerations. The time rate of 
mass flow u in grams per second may be defined by the 
expression 


where p is the density of the gas in grams per cubic 
centimeter, A is the cross-sectional area of the nozzle in 
square centimeters, and v is the velocity in centimeters 
per second. The ideal gas equation may be written in 
the form 


’ 


where P is the absolute pressure in atmospheres, p 
is the density of the gas in grams per cubic centimeter, 
M is the molecular weight of the gas, R is the gas con- 
stant, and 7 is the absolute temperature in degrees 


Kelvin. Equations [8] and [9] may be combined to 
give the expression 
“RT 
A= [10] 


for the area associated with a given mass flow rate at 
any point in the nozzle. The corresponding cross- 
sectional radius is given by the expression 


_ 


From the equation 


which defines the spherical section of the nozzle de- 
picted in Fig. 1, and the fact that y = r, the following 
expression is obtained for the distance downstream 
from the nozzle entrance through the spherical section 


— 9). [13 | 
An expression for the distance downstream through 


the rounded throat section may be derived from the 
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equation which defines the shape of the throat. In 
terms of an 2, y co-ordinate system having its origin 
at the center of the nozzle entrance, the equation of such 
a defining circle is simply 


(x — x)? + (y — yr)? = [14] 


where x, and y; are the co-ordinates of the center of the 
cirele whose radius is 2r, or 2'“r,. Reference to F ig. 1 
shows that 


and that 


= 2 


where 2; is the abscissa of the point at station 1 where 
the circles defined by Equations [12] and [14] are tan- 
gent. By combining Equations [14] and [15] and re- 
membering that y = r, the following expression for the 
distance downstream through the throat section is ob- 
tained: 


z= [2r-2 — (r — [17] 


In Equation [17] the minus sign is valid in the conver- 
gent section of the nozzle and the positive sign in the 
divergent. 

Inspection of Fig. 1 shows that the expression for the 
distance downstream from station 2 is given by the re- 
lation 


[18] 


where rz is the radius of the throat at station 2 where 
the divergent cone of half angle 6 is tangent to the circle 
which defines the shape of the throat. 

Further inspection of Fig. 1 indicates the following 
particular relations. The angle a, associated with the 
point of tangency at station 1, is given by the equation 


cos a = = 0.87868........... [19] 
so that a = 28.516°. Thusat station 1 

= re CoS a = 0.87868r-.............. [20] 

and 

At the throat 

re = = 0.7071 Ire... [22] 

and 

station 2 

re = (3 — 2 cos 15°), = 0.75529r........... [24] 

and 
+ gin 15°r, = 1.5186re.......... [25] 


Kinetic Equations 


The recombination of hydrogen atoms to form hy- 
drogen molecules may be represented by the triple-body 
collision reaction equation 
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where k, and k, are the forward and backward reaction 
rate constants, respectively, and Z is any third body.* 
If [H] denotes the concentration of H in moles per cubic 
centimeter, then the time rate of change of [H] for a 
process in which the volume V is changing is given by 
the expression 
_ 
2 dt 
The validity of the third term on the right-hand side of 
Equation [27] may be easily verified by a considera- 
tion of the identity [1H] = ny/V, where ny is the 
number of moles of H in the volume V of expanding gas. 
In this study it is more convenient to deal with the 
concept of mole fraction Ny than with that of molar 
concentration [H]. By making use of the identities 


1 [H] dV 


= — kel He] +5 (27) 


and 


Equation [27] may be written in the form 


d 


Nu (PV/RT)... . (30) 
For a mass of m grams of gas the ideal gas law may 


be written in the form 


so that 
di In(PV/T) = di In (mR/M) = —- a [32] 
Hence, Equation [30] becomes 


By using Equations [2] and [3] and the value 1.008 
for the atomic weight of hydrogen, the expression for 
the molecular weight of the gas mixture becomes 


M = (2 = Na) 


Equation [33] thus reduces to the following: 


Since the equilibrium constant K, defined by Equa- 
tion [1] and the rate constants in Equation [26] are re- 
lated by the expression 


Equation [35] may be written in the form 
dNu _ Nu) 


dt (RT)? 


(Nu?P? — Nua,P/Kp).....[37] 


Finally, since the velocity v at any distance x down- 
stream is given by the relation 


3 Compare this treatment with that of Penner (3). 
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the rate of change of Ny with distance is given by the 
expression 


aNu k(2 — Nu) y Na P/K,)....%: [39] 


Basic Data 


The pertinent thermodynamic data used in this 
study were taken from the National Bureau of Stand- 
ards (4). They include the total enthalpy above the 
reference temperature 298.16 K and the absolute en- 
tropy for the chemical species H and Hz and the loga- 
rithms of the equilibrium constants for the reaction 
2H = H, at temperatures from 1600 to 3500 K. The 
heats of formation for H and H,2 at 298.16 K are 52,089 
and 0 cal mol, respectively. 

The pertinent physical constants have the following 
values: the acceleration of gravity g = 980.665 cm 
sec~?; the gas constant R = 82.0567 cm* atm deg™! 
mol-!; the mechanical equivalent of heat J = 4.1840 
abs joules cal—. 

The forward reaction rate constant k, in Equation 
[26] is assumed to have the value 10'* mol-? sec! 
over the entire range of temperatures and pressures 
covered by this investigation. This assumption is 
made (a) because the chemical literature is singularly 
lacking in kinetic data for the hydrogen recombination 
reaction‘ and (b) because the primary purpose of this 
study is to show how a rate equation may be dealt with 
in a propellant gas expansion problem. 


Computational Procedure 


In general, the composition of the expanding gas may 
be determined at any temperature and pressure by the 
simultaneous solution of Equations [1] and [2]. The 
entropy of the system, the expansion of which is as- 
sumed to be isentropic, may be determined from the 
initial chamber conditions by means of Equation [4]. 
Once the value of the entropy has been established, 
Equation [4] may be used to determine the pressure at 
any temperature and composition. The velocity of the 
expanding gas at any point may be computed by 
means of Equation [6]. The chamber velocity may be 
determined from the following considerations. 

If, for a given mass flow rate yu, a plot of area A versus 
temperature 7’ is made for various temperatures be- 
tween 7’, and 7, a smooth curve having a well-defined 
minimum is obtained. This point of minimum area 
corresponds to the throat of the nozzle. The chamber 
velocity may, therefore, be computed from Equation 
[10] when the design condition 


is satisfied. 
Since A, cannot be computed explicitly, a tentative 


4 In this connection, see (5) and (6). 


value for A, must be assumed on the basis of some non- 
zero chamber velocity. Then values of A may be com- 
puted for various values of 7. From a plot of these 
data a throat area corresponding to the assumed cham- 
ber area may be determined. If this process is re- 
peated for several values of A, and if a plot of A./A, 
versus A, is made, then, by interpolation, values of A, 
and A, may be obtained which satisfy Equation [40]. 
This procedure for determining A,, A,, and 2, is per- 
fectly general and is independent of the type of flow 
considered. 

The length of the nozzle may be computed from a 
knowledge of end conditions only by means of Equa- 
tions [13], [17], and [18] for constant composition flow 
and for instantaneous chemical equilibrium flow. The 
kinetic case, however, requires a stepwise procedure 
which involves the integration of Equation [39] as well 
as the use of Equations [13], [17], and [18]. For this 
purpose it is convenient to use temperature as a parame- 
ter and to perform the integration over temperature 
intervals which are small enough to insure a degree of 
accuracy consistent with the data used. The two re- 
gions where particular care must be exercised are (a) 
at the very entrance to the nozzle, where the tempera- 
ture gradient is a maximum, and (b) at the throat, 
where the parameters must be determined precisely in 
order to satisfy the design conditions. 

The integration of Equation [39] may be accom- 
plished by writing it in the form 
aNu/dL 


where L denotes the length of the nozzle, and by sub- 
jecting it to an iterative process which makes use of the 
methods and equations of the numerical calculus for 
equal intervals of the argument 7. A list of such equa- 
tions is given by Milne (7). For the present, the fol- 
lowing four-point integration formula is satisfactory 


Ns = No + (3h/8) (No’ + 3N1’ + 3N2’ + Nj’)... . [42] 


In this equation h is the temperature interval and the 
primed quantities are the first derivatives with respect 
to temperature of the function N at four consecutive 
temperatures. 

In order to carry out the integration process ef- 
fectively, it is expedient to construct tables of dN,/dL 
and L values for a judicious spread of Ny values at the 
various temperatures considered. Starting values of 
Ny may be obtained by reference to those calculated for 
the instantaneous equilibrium case. 

Once a set of values of Ny is obtained for the temper- 
ature range from 7, to some arbitrary 7',, corresponding 
values of dN;,/dL and of L may be obtained from the 
previously constructed tables by interpolation. Values 
of dL/dT may be derived from a plot of L versus 7 or 
by means of an appropriate differentiation formula as 
given in (7). Equation [42] may then be used to cor- 
rect the original rough values of Ny as often as neces- 
sary to achieve convergence. 
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Results and Discussion 


The mole fraction composition of the gas entering 
the nozzle at a temperature of 3500 K and a pressure 
of 20 atm is 0.1235 H and 0.8765 He, and its specific en- 
tropy is 22.918 cal deg-! gm~'. With these data as 
initial conditions, a mass flow rate of 1000 gm sec~', 
the nozzle configuration of Fig. 1, and a 2-to-1 chamber- 
to-throat-area ratio, the adiabatic flow of hydrogen to 
an exhaust pressure of 1 atm was studied, assuming 
first, constant composition flow, then instantaneous 
chemical equilibrium flow, and finally, kinetic chemical 
equilibrium flow. The results are presented graphically 
in Figs. 2 through 6. 


Inspection shows that the shape of the kinetic flow 
nozzle is identical with that of the instantaneous flow 
nozzle. This is due to the fact that, since the magni- 
tude 10'* of the forward rate constant k, is for all prac- 
tical purposes infinite, the chamber parameters are 
numerically equal (as far as significant figures are con- 
cerned) for the kinetic and the instantaneous equilibrium 
flow cases. Since the chamber parameters determine 
the geometry of the nozzle, it follows that the shape of 
the nozzles must be identical. 

As was to be expected, the kinetic results are for the 
most part intermediate between the constant composi- 
tion and the instantaneous equilibrium flow results. 
This fact was invaluable in estimating starting values 
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for the iteration process used to integrate Equation 
[39]. A temperature interval of 50 K was employed 
throughout the integration process. A smaller interval 
might have given smoother results, but the increased 
labor would have been greater than the accuracy of 
either the initial data or the final results would warrant. 

Fig. 6 depicts the nozzle configuration for the three 
types of flow investigated. It is interesting to note 
(a) that instantaneous equilibrium flow requires a 
longer (and larger) nozzle than does constant com- 
position flow, and (b) that for kinetic equilibrium flow 
the nozzle is slightly shorter than for instantaneous 
equilibrium flow. 

It may also be noted that nozzle design affects spe- 
cific impulse in the following manner. The specifica- 
tion of a 2-to-1 chamber-to-throat-area ratio results in 
a chamber specific impulse of 141.4 sec for the constant 
composition flow nozzle and 136.0 sec for both the in- 
stantaneous and the kinetic equilibrium flow nozzles. 
It is obvious from Equation [10] that a zero chamber 
velocity would require an infinite chamber area, a fac- 
tor which is not feasible in practice. Furthermore, a 
zero chamber velocity would yield, at an exhaust pres- 
sure of 1 atm, a specific impulse of 823.2 sec instead of 
835.4 for constant composition flow, and 873.8 sec in- 
stead of 884.3 for instantaneous equilibrium flow. It is 
thus evident that nozzle configuration, specifically 
chamber-to-throat-area ratio, has a definite effect on the 
magnitude of specific impulse or jet velocity. 

In the present example the effect of chemical kinetics, 
iLe., the rate of recombination of hydrogen atoms, on 
specific impulse is not very pronounced, as the figures 
884.3 and 879.6 at an exhaust pressure of 1 atm for in- 
stantaneous and kinetic equilibrium flow, respectively, 
indicate. This is due, of course, to the magnitude of 
the forward reaction rate constant k, and to the nature 
of the chemical species. In other expanding propellant 
gases, in which the rate constants associated with the 
pertinent chemical reactions are numerically smaller 
than the hydrogen rate constant and the chemical spe- 
cies themselves are less active than hydrogen, the gen- 
eral effect would be to displace the kinetic results 
toward the constant composition flow values. In such 
cases, even if several differential equations similar to 
Equation [39] are involved, the method of integration 
indicated in this paper would be valid, but the com- 


putation would best be performed on automatic equip- 
ment. 


Conclusions 


The following conclusions may be drawn from the 
foregoing study: 

1 The effect of chemical kinetics on rocket pro- 
pellant gases expanding through a nozzle has been de- 
termined by the simultaneous solution of three sets of 
equations based on thermodynamic, nozzle configu- 
ration, and chemical reaction rate considerations. 

2 The fundamental design calculational require- 
ment, aside from the mass flow rate, is the specifica- 
tion of a finite chamber-to-throat-area ratio. 

3 Nozzles designed on the basis of chemical kinetics 
for perfect expansion are somewhat smaller than in- 
stantaneous equilibrium flow nozzles but are larger 
than constant composition flow nozzles, the difference 
depending on the magnitude of the reaction rates in- 
volved. 

4 The method of stepwise integration described 
may be applied to propellant systems which involve 
several differential equations with, of course, a cor- 
responding increase in computational labor. 
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Jet Propulsion News 


By C. F. WARNER, Purdue University, Associate Editor 


Missiles 


HE U. 8S. Navy’s Martin Viking high-altitude re- 

search rocket set a new world’s record of 135 miles 
for single-stage rockets on Aug. 7, 1951. The Viking 
reached a top speed of 4100 mph. Previous record for 
high altitude rockets was 114 miles reached by the Ger- 
man V-2. The Viking’s liquid oxygen and ethyl alcohol 
rocket engine burned for 75 sec and in that time it had 
reached an altitude of 25 miles from which point it 
coasted the rest of the way. The maximum altitude 
was reached in 4 min and 23 sec. The Viking, 48 feet 
long weighing 5'/. tons, was the seventh of 10 to be 
built by the Glenn L. Martin Company for the Naval 
Research Laboratory. 

++ + 


THE U. S. Air Force will activate the 1st Pilotless 
Bomber Squadron (Light) at the U.S.A.F. Missile Test 
Center, Cocoa, Fla., on Oct. 1, 1951. 

It will be the first operational unit of its kind in the 
Air Force. The unit will be equipped initially with the 
Matador B-61 pilotless bomber, Fig. 3, now in produc- 
tion by The Glenn L. Martin Company, Baltimore, Md. 


FIG. 1. Viking JUST BEFORE FIRING 


The squadron will be supervised in its initial training 
by an already established training unit, the 6555th 
Guided Missile Wing, a part of the Air Research and 
Development Command. Later, it will be made com- 
bat-ready by extensive unit training under the super- 
vision of the Tactical Air Command. 

++ + 


THE Hycon Manufacturing Company and the 
Oederkerk & Ludwig Company of Pasadena, Calif,, 
have a joint multimillion-dollar Navy contract for the 
manufacture of an air-to-air rocket-type missile de- 
veloped by the personnel of the Navy Ordnance Test 
Station in Pasadena. The missile, designed to be 
carried in quantity by high speed jet planes, is powerful 
enough to destroy any known aircraft with one direct 
hit. 


Aircraft 


A U.S. Navy Skyrocket, Fig. 4, piloted by Douglas 
test pilot Bill Bridgeman, has broken its own and all 
existing altitude records. The exact altitude reached by 
the supersonic plane in a flight over Muroe Dry Lake, 
Calif., on August 15, 1951, has not been revealed, al- 


FIG. 2. Viking JUST AFTER FIRING 
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though the officially recognized record is 72,394 ft, 
established in 1935 by a man-carrying balloon. 

The Navy said this latest flight by Mr. Bridgeman in 
the 45-foot rocket-powered plane successfully concludes 
the first phase of a research and development program 
three years old. 

The three Skyrocket planes used in the tests will be 
turned over to the National Advisory Committee for 
Aeronautics for the next phase of a planned research 
program, Rear Admiral T. F. Combs, chief of the Navy 
Bureau of Aeronautics, said. 

The Navy said that the latest flight was made on a 
clear bright day. The Skyrocket, which is powered by 
an assembly of four rocket motors. was air-launched 
from a modified bomber at an altitude of 35,000 ft. 
Within ten sec the plane had exceeded sonic speed and 
was climbing away from the earth at a speed of about 
1000 mph. Almost immediately, the plane was pulled 
into a steeper angle of climb to save fuel for climb rather 
than speed. 

The turbojet engine had been removed to make availa- 
ble space for additional rocket motor propellants. 
During powered flight a ton of propellants per minute is 
expended. 

The Douglas Aircraft Company designed and built 
the Skyrocket to meet the special requirements for high 
speed and extreme altitude research established by the 
Navy Bureau of Aeronautics and the National Ad- 


FIG. 3. U. 8S. AIR FORCE B-61 MARTIN Matador p1LOTLESS BOMBER 


visory Committee for Aeronautics to enable procure- 
ment of much needed aerodynamic information in the 
transonic and supersonic speed ranges. The plane was 
equipped with four hundred openings on the wing and 
tail surfaces for air pressure determinations, and over 
25 miles of electrical wiring were used to connect the 


electrical instruments. 


MecDONNELL Aircraft Corporation has announced 
that its new Navy fighter, the XF3H-1 Demon jet air- 


craft has been flight tested for a period of 24 minutes. 
++ + 


THE Glenn L. Martin Company has just received an 
English Electric Canberra light bomber from England. 
This plane will be built in substantial quantities by the 
Martin Company for the U. S. Air Force. Kaiser 
Metal Products, Inc., of Bristol, Pa., is a subcontractor 
of the Martin Company in connection with this pro- 
gram. 

Rocket Engines 


THE Curtiss-Wright Corporation’s XLR-25 rocket 
engine, the first rocket power plant capable of being 
throttled from full thrust to 25 per cent full thrust to be 
built in this country, has passed reliability tests at the 
corporation’s Propeller Division Plant in Caldwell, 
N. J. The motor should be ready for acceptance tests 
by the Air Force soon. This motor is designed for use 
in the Bell Aircraft Corporation’s X-2. 


FIG, 4. NAVY’S RESEARCH AIRPLANE D-558-2 Skyrocket 
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Turbojet Engines 


THE Pratt & Whitney Aircraft J-42 Turbo-Wasp has 
become the first jet aircraft engine to reach an officially 
authorized service life of 1000 hours between major over- 
hauls, it was announced recently by Willam P. Gwinn, 
general manager. The new 1000-hour overhaul interval 
for the J-42 brings this Pratt & Whitney jet engine to a 
level of operational reliability and durability compara- 
ble to piston engines. The J-42 is a centrifugal flow 
type jet engine developed by Pratt & Whitney Aircraft, 
under auspices of the Navy’s Bureau of Aeronautics, 
from an original design of Rolls-Royce, Ltd. It is 
rated at 5000 Ib static thrust dry and 5750 lb using 
water injection. 


New Facilities 


THE U. 8. Navy’s Bureau of Yards and Docks is 
completing an aeronautical turbine laboratory at 
Trenton, N. J. Incorporating designs of the York 
Corporation, the laboratory will consist of five testing 
areas composed of two altitude “strato-chambers,” two 
test cells, and one turboprop cell. The world’s largest 
liquid rheostats, to help test full size jet engines, guided 
missiles, wing sections, and fuselages, are to be built 
by the Westinghouse Electric Corporation. 


+ 


PRESIDENT Harry Truman recently dedicated the 
Arnold Engineering Development Center at Tullahoma, 
Tenn. The $1,575,000,000 institution, to be operated 
by a private corporation, Aro, Inc., will begin test 
operations early next year. The facilities will include 
turbojet and ramjet test units, gas dynamics facilities 
for testing models at speeds above the Mach 5 range, 
and propulsion wind tunnels with transonic and super- 
sonic test sections. 

++ + 


THE recently established National Bureau of Stand- 
ards laboratory at Corona, Calif., will be devoted pri- 
marily to Defense Department research projects. 
R. D. Huntoon, who has been named associate director 
to head the new laboratories, said that the most 
important activity there will be the development of 
guided missiles. Every phase of guided missile develop- 
ment from basic and applied research to the fabrication 
of experimental parts will be carried on. An analogue 
computer which will compute trajectories of proposed 
missiles will be installed. Other facilities at the former 
Navy installation will include electronic laboratories, 
machine shops, a wind tunnel, jet engine test cells, 
altitude chamber, a missile assembly section, and a 
technical library. 


+ 
LOCKHEED Aircraft Corporation reported that it 


will construct a $12,615,000 aircraft assembly plant for 
the U.S. Air Force at Palmdale Airport, 65 miles north- 
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east of Los Angeles, Calif. The plant will employ several 
thousand persons in the assembly of jet aircraft. 


Air Research and Development 
Command 


FUTURE plans for guided missiles, supersonic planes, 
and other technological advancements are being pre- 
pared at the Baltimore headquarters of the new Air 
Research and Development Command. The new unit, 
known as ARDC for brevity, is dedicated to maintaining 
continuous United States air supremacy by exploring 
the unlimited boundaries of scientific knowledge. 
ARDC headquarters will direct operations at proving 
grounds throughout the nation. ARDC will eventually 
have assigned to it all research and development activi- 
ties of the Air Force. 

The new command is the outgrowth of the findings 
of a committee headed by Dr. Louis N. Ridenour, the 
dean of the University of Illinois graduate school. In 
1949, General Hoyt 8. Vandenberg, Chief of Staff, re- 
quested that such a committee conduct an intensive 
survey of Air Force research and development organiza- 
tion, procedures, policies, and facilities. The Ridenour 
committee recommended that a new command be 
formed to be responsible for all Air Force research and 
development. 

It was announced on May 23 that Lt. General 
Earle E. Partridge is to succeed General Schlatter as 
ARDC commanding general. With the initial phase 
of the new command’s development nearing its close, 
General Schlatter has been assigned as Commanding 
General of the Allied Air Forces in Southern Europe. 
Brigadier General Donald J. Keirn is director of re- 
search, and Brigadier General John W. Sessums, Jr., is 
director of operations. Dr. Mervin J. Kelly and 
Donald A. Quarles, top-ranking civilian scientists are 
serving as part-time advisers to General Schlatter and 
will continue in that capacity. 


Personalities 


COL. KENNETH E. FIELDS has been appointed 
director of the Division of Military Application, U. S. 
Atomic Energy Commission, succeeding Brig. Gen. 
James McCormack, Jr., who has been assigned as a 
special assistant to deputy chief of staff for development 
in U.S. Air Force headquarters. 

+~+ + 

REAR ADM. THOMAS SELBY COMBS, U.S.N., 
took over the duties of chief of the Bureau of Aero- 
nautics, Department of the Navy, on May 1, 1951, fora 
four-year term. Prior to this assignment, he was chief 
of staff and aide to commander-in-chief, Atlantic Fleet. 
He had been with Bureau of Aeronautics previously, 
serving from 1946 to 1948 as its deputy chief. 

+ 

D. R. SHOULTS, currently the director of engineer- 
ing for Aro, Inc., St. Louis, Mo., has been appointed 
director of the General Electric Company’s Aircraft 
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Nuclear Propulsion Project for the Air Force and 
Atomic Energy Commission, C. W. Lapierre, manager 
of the company’s Aircraft Gas Turbine Division, an- 
nounced recently. The project is engaged in the fur- 
ther development of an atomic-powered engine for air- 
planes. Headquarters of the project will be at the 
company’s Lockland, Ohio, turbojet research, engineer- 
ing, and manufacturing center. 
++ + 

ELMER E. NELSON moved to his new quarters in 
Sacramento, Calif., to assume his new duties as resident 
manager of the Aerojet Sacramento Division. The new 
rocket plant, located 16 miles northeast of Sacramento, 
was constructed for the purpose of augmenting the 
present Azusa, Calif., facilities in the production of 
solid propellant rocket motors. 

++ + 

C. D. PERRINE, JR., former specialist in missile 
homing guidance systems, has been named assistant 
chief engineer in charge of electronics at Convair, San 
Diego. 
++ + 


JOINT announcement has been made by Harry 
F. Guggenheim, president of the Daniel and Florence 
Guggenheim Foundation, Harold W. Dodds, president, 
Princeton University, and Lee A. DuBridge, president, 
California Institute of Technology, of the appointment 
of seven new Daniel and Florence Guggenheim Jet 
Propulsion fellows. The fellows, selected on the basis of 
their technical ability, deep interest in rockets and jet 
propulsion, and demonstrated leadership qualities, will 
carry out programs of advanced study and pioneering 
research at the Daniel and Florence Guggenheim Jet 
Propulsion Centers established by the Foundation at 
Princeton and Caltech late in 1948. 

Recipients of the grants for study at Princeton are: 
Arnold Goldburg, Mahwah, N. J., a graduate student at 
Massachusetts Institute of Technology; Marcus P. 
Knowlton, II, Auburndale, Mass., now studying at 
Princeton University; Lt. D. M. Saunders, Takoma 
Park, Md., at present attending the U. S. Naval Post- 
graduate School at Annapolis, Md.; and Marcel 
Vinokur, Long Island City, N. Y., a senior at Cornell 
University. 

The fellowships at the Caltech Center are awarded 
to: George E. Hlavaka, Racine, Wis., now studying 
toward a Ph.D. degree at Caltech; Robert C. Evans, 
Pittsburgh, Pa., at present a student in jet propulsion 
at Caltech; and Herman Carl Thorman, Berea, Ohio, a 
graduate student in the School of Aeronautics, Purdue 
University. 

The fellowships carry a stipend of $2000 annually, 
and are normally awarded for one year, but may be 
renewed for an additional year. All the Princeton 
fellows appointed last year have had their grants re- 
newed. These are Major Dale D. Davis, U.S.A.F., 
Effingham, IIl.; Clarence B. Cohen, Mt. Vernon, N. Y.; 
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John E. Scott, Jr., Portsmouth, Va.; and David H. 


Ross, Chicago, Ill. Caltech fellows who have had 
their grants renewed are: Thomas C. Adamson, Jr., 
Lagrange, IIl.; Eldon L. Knuth, Launa, Iowa; and 
Edward E. Zukoski, Birmingham, Ala. 
+ 

THE Rocket Physics Department of the Aerojet 
Engineering Corp., as of July 18, 1951, is no longer a 
separate department. A consolidation has been effected 
between this department of the Research Division and 
the Liquid Engine and Proving Grounds Departments. 
The over-all supervision of the Liquid Engine Depart- 
ment remains in the hands of C. C. Ross, chief engineer, 
with R. B. Young as assistant chief engineer. 

+ 

RECENT changes in the executive staff of Reac- 
tion Motors, Inc., were announced by Charles W. New- 
hall, Jr., Executive Vice President. R.W. Young, for- 
merly Vice President in charge of Engineering at 
Wright Aeronautical Corporation, has been elected 
President of R.M.I., replacing Lovell Lawrence, Jr., 
who resigned. H. B. Horne, Jr., has been appointed 
Manager of Engineering, and W. P. Munger, Chief 
Engineer, replacing H. R. Moles who resigned. 


+ 


ERNEST J. VOGT 


One of the purposeful young men of rocketry, Ernest J. 
Vogt, passed on to his Creator during the early hours of 
September 27, 1951, at the age of forty. Although ill for 
many months, his mental faculties were unimpaired until 
he died, and to his great and constant satisfaction he was 
able to work on special development projects until the end. 

Mr. Vogt began his technical career as a communications 
engineer, and had already achieved wide recognition for 
his contributions in this field when he was induced to join 
the starting nucleus of the Aerojet Engineering Corpora- 
tion in 1942, thus launching a new career in rocket develop- 
ment. He was quickly promoted to Assistant Chief En- 
gineer, but he later found that the most satisfactory means 
of expressing his high creative abilities was in the field of 
heading up self-contained special projects. Of the many 
important projects that he carried out, there is one that is 
sufficient in itself to earn him a permanent place on the 
honor roll of rocketry. He was the Aerojet engineer in 
charge of designing, building, and finally flight testing the 
power plant for the XP-79, the first United States airplane 
powered solely by rocket motors and actually flown. This 
was in the fall of 1943. 

His interest in the broad aspects of rocket propulsion 
was always keen and imaginative, and he was a constant 
source of stimulation and enthusiasm to those who were 
fortunate to be his colleagues. It is reported by Mr. 
Andrew G. Haley, executor of Mr. Vogt’s estate, that 
among the personal effects by his bedside when he died was 
an unfinished manuscript on the technical problems of 
space travel. 

He leaves his two children, Andrew John, eight years 
of age, and Mary Michaela, aged six, and a great many 
personal friends in the American Rocket Society and else- 
where who mourn his passing. 
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American Rocket Society News 


By A. F. BOCHENEK, Associate Editor 


Record American Rocket Society 
Program for 1951 National Convention 


vas record program of the 1951 Na- 
tional Convention to be held at the 
Chalfonte-Haddon Hall Hotel, Atlantic 
City, N. J., November 28-30, 1951, is 
ample evidence that the American Rocket 
Society provides an acknowledged na- 
tional forum for men who are contributing 
to the field of jet propulsion and rocket 
technology. 

The three-day program—consisting of 
five technical sessions at which 17 papers 
will be presented, an annual meeting 
of members, and a social program— 
is expected to attract many of the key 
figures in the rocket industry and in jet 
propulsion. 

Because the ARS is an affiliate of The 
American Society of Mechanical Engi- 
neers, ARS members will be able to attend 
without paying any registration fee, the 
1951 annual meeting of the ASME which 
is being held concurrently in the same 
hotel. The diversified ASME program 
touching on many engineering problems 
offers much of common interest to ARS 
members. 


Annual Meeting 


The ARS annual meeting will be held on 
Wednesday, at 9:30 a.m. President H. R. 
J. Grosch will open the meeting with a re- 
port of the Board of Directors. He will be 
followed by G. Edward Pendray, ARS 
treasurer, who will report on Society 
finances. A progress report on the year’s 
activities in membership development 
will be presented by R. W. Porter, chair- 
man of the Membership Committee. 

Following election of a president, vice- 
president, directors and nominating com- 
mittee for 1952, the meeting will hear an 
account of the Second International Con- 
gress on Astronautics held recently in 
London, England, from Lt. Comdr. F. C. 
Durant III, who was an official ARS 
delegate to the Congress. 


ARS Section Luncheon 


At the ARS Section Luncheon at noon 
Wednesday, important administrative and 
organizational problems will come under 
discussion. The luncheon is open to Sec- 
tion officers and official delegates of the 
individual sections. President Grosch will 
preside. He will report briefly on 1951 
activities and will introduce official rep- 
resentatives of the following ARS Sec- 
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tions, who will report on local activities: 
New York, Indiana, Southern California, 
and New Mexico-West Texas. 

The rapid growth of the rocket industry 
and the large number of engineers that are 
being enlisted to attack various rocket 
problems have created a need for more 
ARS Sections to serve the technical needs 
of these widely scattered centers of rocket 
interest. This problem will be discussed 
by R. W. Porter, chairman of the Member- 
ship Committee. The advisability of as- 
signing the entire membership on a geo- 
graphical basis to Sections already organ- 
ized will be introduced for comment from 
the floor. 

Of immediate interest to all members 
will be the comments of Martin Summer- 
field, Editor-in-Chief, on development of 
the new ARS Journal. 

The luncheon program will end with a 
general discussion of future activities, es- 
pecially the advisability of planning 
similar luncheons at future national con- 
ventions of the Society. 


Honors Night 


The high light of the convention will be 
the banquet and honors night on Thurs- 
day when the Society will pay tribute to 
creative engineers and scientists who by 
their work have advanced the technology 
of rocket and jet propulsion. 

The main speaker will be the Honorable 
A. 8. Alexander, Under-Secretary of the 
Army, who will discuss “Matters Relating 
to National Defense.” 

In a special ceremony awards will be 
presented as follows: 


1. Goddard Memoria! Lecture Award 
to Comdr. Robert C. Truax, U.S. 
Navy. 

2. C. N. Hickman Award to Dr. 
William H. Avery, Supervisor, 
Propulsion Group, Applied Phys- 
ics Laboratory, Johns Hopkins 
University, Silver Spring, Md. 

3. G. Edward Pendray Award to 
George P. Sutton, consultant, 
North American Aviation, Los 
Angeles, Calif. 

4. ARS Student Award to David Elliott 
and Leo Rosenthal. 

Fellow memberships in the Society, ac- 

knowledging eminence in the field of 


rocket and jet propulsion engineering, will 
be bestowed on: 

Chandler C. Ross, Aerojet Engineering 
Corporation, Azusa, Calif. 

Benjamin F. Coffman, Jr., Ship Insti|- 
lations Division, Bureau of Aero- 
nautics, Navy Dept., Washington, 
D. C. 

E. H. Hull, General Electric Company, 
Schenectady, N. Y. 

Charles E. Bartley, Jet Propulsion 
Laboratory, California Institute of 
Technology, Pasadena, Calif. 

Col. H. N. Toftoy, chief, Rocket Branch, 
Research and Development Division, 
Office of Chief of Ordnance, Washing- 
ton, D. C. 


Technical Program 


The tentative program foliows: 


3:00 p.m., Wednesday, November 28 

Chairman: C. C. Ross, chief engineer, 
liquid engine department, Aerojet Engi- 
neering Corporation, Azusa, Calif. 

Vice-Chairman: Powel Brown, rocket 
test engineer, special projects depart- 
ment, M. W. Kellogg Company, Jersey 
City, N. J. 

Attitude Stabilization for Supersonic Ve- 
hicles, by C. W. Besserer, supervisor of 
engineering, and A. J. Bell, staff engi- 
neer, Applied Physics Laboratory, The 
Johns Hopkins University, Silver Spring, 
Md. 

Some Experimental Dynamic Launching 
Techniques for Testing Aircraft 
Rockets, by A. W. Nelson, Nava) Ord- 
nance Test Station, Inyokern, Calif. 

Unusual Applications of the Momentum 
Principle, by T. F. Reinhardt, U. S. 
Naval Air Rocket Test Station, Lake 
Denmark, Dover, N. J. 


9:30 a.m. Thursday, November 29 

Chairmen: J. H. Sheets, chief rocket de- 
sign engineer, Curtiss-Wright Corpora- 
tion, Propeller Division, Caldwell, N.J. 

Vice-Chairman: M. Meyer, supervising 
project engineer, Curtiss-Wright Cor- 
poration, Propeller Division, Caldwell, 
N. J. 

Installation of Rocket Engines in Air- 
planes, by F. A. Coss, Reaction Motors, 
Inc., Rockaway, N. J. 

Description of the NOTS Aeroballistics 
Laboratory, by I. E. Highberg, Naval 
Ordnance Test Station, Inyokern, Calif. 

The United States Air Force Experi- 
mental Rocket Engine Test Station, 
Edwards, Calif., by Lt. Richard A. 
Schmidt, chief, test engineering section, 
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Rocket Branch, and Donald L. Dynes, 
rocket engine test engineer, Rocket 
Branch, Edwards Air Force Base, 
Muroc, Calif. 


2:30 p.m. Thursday, November 29 

Chairman: Col. N. C. Appold, U. S. 
Air Force, chief, Power Plant Labora- 
tory, Aeronautics Division, Wright Air 
Development Center, Dayton, Ohio. 

Vice-Chairman: Lt. Col. E. N. Hall, U.S. 
Air Force, assistant chief, nonrotating 
engine branch, Power Plant Laboratory, 
Wright Air Development Center, Day- 
ton, Ohio. 

Isothermal Combustion Under Flow Con- 
ditions, by Capt. J. A. Bierlein and K. 
Scheller, power plant research engi- 
neers, Wright Air Development Center, 
Wright-Patterson Air Force Base, Day- 
ton, Ohio. 

Porous, Sweat and Film Cooling With 
Reactable Coolants, by Luigi Crocco, 
Robert H. Goddard Professor, The 
Daniel and Florence Guggenheim Jet 
Propulsion Center, Princeton University, 
Princeton, N. J. 

High Flux Heat Transfer and Coke 
Deposition of JP-3, by J. B. Hatcher 
and D. R. Dartz, Jet Propulsion Labor- 
atory, California Institute of Technol- 
ogy, Pasadena, Calif. 

The Principle of the Concentric Nozzle, by 
W. F. Kaufman and B. N. Abramson, 
U. S. Naval Air Rocket Test Station, 
Lake Denmark, Dover, N. J. 


9:30 a.m. Friday, November 30 

Chairman: C. M. Hudson, technical as- 
sistant to the chief, Rocket Branch, Re- 
search & Development Division, Office 
of Chief of Ordnance, Washington, 
D.C. 

Vice-Chairman: H. F. Calcote, Experi- 
ment, Inc., Richmond, Va. 

Combustion Studies with a Rocket Motor 
Having a Full Length Observation 
Window, by Kurt Berman and Stanley 
Kk. Logan, General Electric Company, 
Malta Test Station, Ballston Spa, 
N.Y. 

Photographic Techniques Applied to 
Combustion Studies: — Two-Dimen- 
sional Transparent Thrust Chamber, by 
J. D. Thackrey and John H. Altseimer, 
development engineer, Aerojet Engineer- 


ing Corporation, Azusa, Calif. 

Experimental Problems in High Pressure 
Combustion, by R. L. Wehrli, chief 
physicist, Reaction Motors, Inc., Rocka- 
way, N. J. 


2:30 p.m. Friday, November 30 
Chairman: George P. Sutton, supervisor, 
Aerophysics Laboratory, North Amer- 
ican Aviation, Inc., Downey, Calif. 
Vice-Chairman: O. K. Doyle, specialist, 
Aerophysics Laboratory, North Ameri- 
can Aviation, Inc., Downey, Calif. 

low Stability in Small Orifices by R. P. 
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Northrup, General Electric Company, 
Malta Test Station, Ballston Spa, 
N: 


Injector Spray and Hydraulic Methods of 


Rocket Motor Analysis, by Kurt R. 
Stehling, rocket research engineer, Bell 
Aircraft Corporation, Buffalo, N. Y. 
Fluctuations in a Spray Formed by Two 
Impinging Jets, by Marcus F. Heid- 


mann, aeronautical research scientist, 
and Jack C. Humphrey, NACA Flight 
Propulsion Laboratory, Cleveland, 
Ohio. 

Some Interactions of Circular Liquid Jets, 
by K. Dexter Miller, Jr., project engi- 
neer, special projects department, The 
M. W. Kellogg Company, Jersey City, 
N. J. 


R. W. Porter Reviews V-2 


Launching Program at 


New York Chapter Meeting 


The work horse of high altitude re- 
search—the German-designed V-2—came 
in for some well-earned plaudits recently 
from Dr. R. W. Porter who is in charge of 
the U. S. Army’s Guided Missile Project 
for the General Electric Company. With 
vehicles capable of carrying payloads as 
large as 4000 lb, the V-2 program made it 
possible for the personnel of two rocket 
proving grounds, 15 university research 
groups, and many government and indus- 
trial designers to “cut their teeth” in the 
relatively new art of building and launch- 
ing rockets. . 

Dr. Porter spoke before 100 ARS mem- 
bers and guests of the New York Chapter 
of the American Rocket Society at its 
monthly meeting on October 19, 1951. 
While his topic was ““V-2 Troubles’ he 
made it abundantly clear that even from 
the launching failures valuable lessons 
were learned. 


With the aid of German engineers and 
technicians, who eagerly turned over to 
American engineers available parts and 
records of the German V-2 program, the 
Army was able to launch successfully 45 
V-2’s out of 68 firings. Of these, 32 were 
perfect launchings with no malfunctioning 
of any of the components. The remaining 
13 showed minor erratic behavior after 
take-off, although not serious enough to 
interfere with the research purposes of the 
firings. 

This record, Dr. Porter said, was con- 
siderably better than that originally pre- 
dicted by the German eng neers. Because 
of better workmanship in some of the 
American-produced components, the V-2’s 
reached higher altitudes and carried 
heavier payloads than the V-2’s fired in 
Germany, although the units were of 
exactly the same design. 

In justifying the cost of the program, Dr. 


(Left to right: L. D. White, director of Hermes Project, White Sands Proving Grounds, 


and director of ARS New Mexico Section. 


C. W. Chillson, consulting engineer, Curtiss 


Wright Corporation, and vice-president of ARS. R. W. Porter, division engineer, General 
Electric Company, Schenectady, N. Y., and director of ARS. Lee J. Bregenzer, design 
engineer, Port of New York Authority, New York, N. Y., and director of ARS New 
York Section. Roy Healy, division head, operating division, M. W. Kellogg Company, 
Jersey City, N. J., and past president of ARS.) 
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Porter pointed out that the firings saved 
the nation time and money in training re- 
search personnel, firing crews, and in- 
dustrial designers, and produced some in- 
teresting and useful upper atmosphere 
data which could be obtained in no other 
way than by full-scale tests on actual 
rockets of the size of the V-2. 

Dr. Porter praised the aid given the 
Government by the 15 American univer- 
sities who voluntarily organized an Upper 
Atmosphere Research Panel to advise on 
the use of the V-2’s. Government agen- 
cies gladly accepted the aid of the pane! 
and turned over to it the responsibility of 
allocating payloads of successive V-2’s 
to various research projects. 

Eight of the firings were of the two-stage 
type, using a V-2 to launch a smaller Wac 
Corporal. The perfection of a multiple 


staging technique is considered one of the 
most significant results of the program. 
But the record of achievement, made pos- 
sible by the large V-2 payloads and the 
ingenuity of the co-operating research 
teams, contained a long list of aerody- 
namic, heat-transfer, and component- 
design data. 

In the rundown of the firing failures, 
Dr. Porter attributed 22 to poor design, 
wearout, and poor workmanship of com- 
ponents, 8 to inexperience in launching 
technique, 4 to special experiments con- 
ducted in connection with the firings, 1 to 
personnel laxity, and 1 to a random cause. 
Dr. Porter expressed confidence that in a 
like number of future firings no more than 
5 failures would occur. 

The evening ended with a short film on 
V-2 operations and a discussion period 
followed by refreshments. 


ARS Joins International Astronautical 
Federation 
By Lt. Comdr. F. C. Durant III (U.S.N.) 


N September 3-8, 1951, the Second 
International Astronautical Congress 
was held in London, England. Hcst for 
the week-long convention was the British 
Interplanetary Society. The purpose of 
this Congress was to establish an Inter- 
national Astronautical Federation and to 
present technical papers on the problems of 
realizing earth satellite vehicles and inter- 
planetary travel. Presiding at the Con- 
gress was Provisional President Dr. 
Eugen Sanger who was elected as first 
President of the International Astronau- 
tical Federation. Dr. G. Loeser (G.f.W.) 
and Andrew G. Haley (ARS) were elected 
Vice-Presidents for matters relating to 
Europe and the Unite ’ States, respectively. 
Also attending the Congress were Prof. 
Herman Oberth and Dr. Irene Bredt. 
In addition to the United States and 


Great Britain, authorized delegates were 
present from societies in France, Sweden, 
Germany (British and U. S. Zones), Aus- 
tria, Switzerland, Italy, Spain, and Argen- 
tina. Representing the American Rocket 
Society, Andrew G. Haley, as General 
Counsel, was principal delegate on par- 
liamentary matters, and Lt. Comdr. Fred- 
erick C. Durant III (U.S.N.) acted as 
principal delegate on technica] matters. 

Lt. Comdr. Durant presented a paper 
for ARS member Dr. Wernher von 
Braun, entitled “The Importance of 
Satellite Vehicles in Interplanetary Flight.” 
The subject of this paper was the keynote 
of the Congress. The importance of 
the artificial satellite as a take-off point 
for interplanetary travel received com- 
plete acceptance by the assembled dele- 
gates. Dr. von Braun’s paper was actu- 


PRINCIPAL DELEGATES PRESENT AT SECOND INTERNATIONAL ASTRONAUTICAL CONGRESS 


Country Scciety 
Great Britain British Interplanetary 
Society 
Germany Gesellschaft fiir 
Weltraumforschung 
Nordwestdeutsche G.f.W. 
Austria Osterreichische G.f.W. 
France Groupement Astronautique 
Frangais 
Sweden S.LS. 
Switzerland A.S. fiir Weltraumfor- 
schung 
Italy A.LR. 
Spain A.E.A. 
Argentina Sociedad Argentina Inter- 
planetaria 
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Location Principal Delegates 

London Arthur C. Clarke 

A. V. Cleaver 
Stuttgart Dr. G. Loeser 

Heinz H. Kdlle 
Frankfurt Albert Stemmer 

F, K. Jungelaas 
Innsbruck Friedrich Hecht 
Paris A. Ananoff 
Stockholm Ing. Ake Hjerstrand 
Baden Ing. Josef A. Stemmer 
Bologna Dr. Giovanni Casiraghi 

Prof. Eula 

Capt. F. Partel 
Madrid Tomas Mur 


Buenos Aires Prof. Teofilo M. Tabanera 


ally a synopsis of a book, entitled “Mars 
Project,’ to be published shortly. The 
scientific appendix of this work presents 
data and calculations showing that the 
artificial satellite and subsequent cir- 
cumplanetary flight are technically feasible 
today if money, personnel, material, and 
technical direction were made available. 
This realization that the manned, arti- 
ficial satellite is a logical and expected step 
beyond current rocket-powered missi!es 
has never before been so universally ac- 
cepted by professional engineers and scien- 
tists. 

The Detroit Rocket Society, Pacilic 
Rocket Society, and Reaction Research 
Society were represented by E. V. Sawyer 
and David Elliott. 

Steps were taken to standardize nomen- 
clature and symbols relating to rocket 
propulsion and astronautics. In addition, 
a home address in Switzerland is being ob- 
tained. An international archive has been 
established to provide reference material 
on subjects relating to astronautics. Ing. 
Josef Stemmer of Baden will coordinate 
this effort. Vice-Presidents Loeser and 
Haley were directed to prepare a 
draft constitution, which will be circu- 
lated to the member societies prior to the 
next meeting. The next meeting of the 
Federation will be in 1952 at Stuttgart. 


@ MICROTORQUE Variable Resistors and 
Potentiometers require as little as .003 
in. oz. torque to operate. This unique fea- 
ture makes the MICROTORQUE invaluable 
for applications where the position of in- 
strument pointers, gyroscopes, and delicate 
instruments in general must be recorded, 
transmitted or indicated at a distance, 
and Giannini are the sole makers of 
MICROTORQUE Potentiometers. 

A variety of resistance values and 
circuits available. 


Write for booklet. 
G.M. Giannini & Co., Inc. 
Pasadena 1, California 
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Technical Literature Digest 


By H. S. SEIFERT, California Institute of Technology, Associate Editor 


CoNTRIBUTORS: 


D. Altman 
D. I. Baker 
R. B. Canright 


F. C. Gunther 
S. S. Penner 
H. S. Seifert 


Eprror’s Note: 

The following collection of references is 
not intended to be comprehensive, but is 
rather a selection of the most significant 
and stimulating papers which have come 
to the attention of the contributors. The 
reader will understand that a considerable 
body of literature is unavailable because 
oi security restrictions. We invite con- 
tributions to this department of references 
which have not come to our attention, as 
well as comment on how the department 
may better serve its function of providing 
leads to the jet propulsion applications of 
many diverse fields of knowledge. 


Book Reviews 


Rockets, Missitks AND Space TRAVEL, 
by Willy Ley, The Viking Press, New 
York, 1951, xii + 436 pp. $5.95. 


Reviewed by H. 8. Serrert 


HE popular appeal of the subject 
of rockets and space travel is demon- 
strated by the fact that this volume of 
Mr. Ley’s is essentially a sixth edition of 
“Rockets,” first published in 1944. Span- 
ish, Italian, and German editions also 
exist. 

“Rockets, Missiles and Space Travel” 
is a history of rocketry, compiled in 
roughly chronological order, which de- 
scribes the efforts and hardships of the 
early visionaries and rocket experimental- 
ists, and which shows how their dreams 
have been in part fulfilled. It is written 
with great sympathy and insight into the 
lives and thoughts of these men, and is 
well garnished with incidents calculated 
to emphasize the human and dramatic as- 
pects of their lives. Mr. Ley is in an ad- 
vantageous position to do this since he 
himself was a charter member and officer 
of the original German rocket society, the 
VfR (Verein fiir Raumschiffahrt). 

The book contains chapters on the early 
classical period of the concept of space 
flight, the Jules Verne or fantasy era, the 
first weapon rockets, and the Russian, 
French, German, and American pioneers. 
This material is followed by three chap- 
fers new to this edition: One on weapon 
tockets, one on the German rocket de- 
velopment center at Peenemunde, and the 
last on the corresponding American center 
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at White Sands Proving Ground in New 
Mexico. The concluding chapters on 
orbital rockets, the space ship, and the 
space terminal are largely discussions of 
the dynamics of long-range flight, with 
little emphasis on the speculative subjects 
of vehicle design or pilot physiology. 

The 90 pages of appendixes contain an 
interesting assortment of relatively new 
information, including a discussion of 
rocket airplanes, the characteristics of 
various missiles, and tables relating to 
space travel such as meteorite and other 
astronomical data. There is a reasona- 
bly complete bibliography of rocket liter- 
ature, including a section on imaginative 
literature! 

The book is popular rather than tech- 
nical in content, written with emphasis 
on the dramatic, and in a discursive and 
at times even rambling style. Neverthe- 
less, it will interest sericus workers in 
rocketry for the wealth of incidental in- 
formation and anecdote which it supplies. 
Physical concepts are accurately and 
soundly discussed, and illustrations are 
plentiful. The reviewer is impressed by 
the author’s evident enthusiasm and the 
breadth of his information. While this 
volume cannot be considered to be a text- 
book of rocket principles, it nevertheless 
should be available in every jet propul- 
sion library as a historical source book. 


THE STRUCTURE AND MECHANICAL PROP- 
ERTIES OF Metats, by Bruce Chalmers, 
John Wiley & Sons, Inc., New York, 
N. Y., 1951, ix + 132 pp. $3.50. 


Reviewed by LAwrENcE J. Hutt 


THE object of the author in “The 
Structure and Mechanical Properties of 
Metals” is to provide a simple, easily 
understood picture of the structure of 
metals and the relation of the structure 
to the mechanical properties. Dr. Chal- 
mers has accomplished this objective to a 
high degree without resort to mathe- 
matical considerations. He explains the 
theory of structure of pure metals, then 
enlarges the picture to include alloys. 
This is supplemented by chapters describ- 
ing the effects of distortion and heat 
treatment on structure. A short chapter 
is included on the general determination 
of structure, followed by an excellent 
treatment of mechanical properties and 


their dependence on structure. This 
L 


chapter should be quite useful for the 
engineer who desires an easily under- 
stood and well presented source of infor- 
mation on what mechanical properties are 
and what they mean. 

The book adheres to the material de- 
scribed by the title and is not a work on 
general metallurgy nor a reference hand- 
book. Inaccuracies may be noted by the 
careful reader, some of which appear in 
the tabular data. It is assumed that the 
reader has a general background in phys- 
ics and chemistry and that many com- 
mon engineering terms will be understood. 
The book is recommended to the engi- 
neer who does not have a general know]- 
edge of the subject and who is seeking a 
readable account of it. It is recom- 
mended to the practicing metallurgical or 
mechanical engineer only for review. 


Books 


Introduction to the Transfer of Heat 
and Mass, by E. R. G. Eckert, MceGraw- 
Hill Book Co., Inc., New York, N. Y., 
1950. 

The Handbook of Measurement and 
Control, by M. F. Behar, The Instruments 
Publishing Co., Inc., Pittsburgh, Pa., 
1951. 

Flow Measurement with Orifice Meters, 
by R. F. Stearns and others, D. Van Nos- 
trand Co., Inc., New York, N. Y., 1951. 

Guided Missiles: Rockets and Torpe- 
does, by Frank Ross, Lothrop, Lee and 
Shepard Co., Inc., New York, N. Y., 
1951, 186 pp. 

The Physics of Powder Metallurgy (a 
symposium), by Walter E. Kingston, 
McGraw-Hill Book Co., Inc., New York, 
N. Y., 1951, 404 pp., $8.50. 

Handbook of Supersonic Aerodynamics, 
NAVORD Report 1488, The Johns Hop- 
kins Univ., APL, Silver Spring, Md., 
July 15, 1951. 


Jet Propulsion Engines 


Investigations of Axial Flow Com- 
pressors, by J. T. Bowen, R. H. Sabersky, 
and W. D. Rannie, Transactions of The 
American Society of Mechanical Engineers, 
vol. 73, January, 1951, pp. 1-15. 

U.S. Research Project SQUID: Rocket, 
Ramjet and Pulse-Jet, by Hugh Harvey, 
Shell Aviation News, August, 1950, pp. 
16-18. 
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Burners for Supersonic Ramijets, by 
J. B. Fenn, H. B. Forney, and R. C. 
Garmon, Industrial and Engineering Chem- 
istry, vol. 43, July, 1951, pp. 1663-1672. 

The Solid Fuel Ramjet, by E. F. 
Chandler, Aero Digest, vol. 61, September 
1950, pp. 19, 86, 87. 

Air Inlets, by Robert McLarren, Aero 
Digest, vol. 62, April, 1951, pp. 75-85. 

An Analysis of Helicopter Propulsion 
Systems, by J. A. Johnson, Transactions of 
The American Society of Mechanical 
Engineers, vol. 73, July, 1951, p. 519. 

The Generalized Approach to the Se- 
lection of Propulsion Systems for Air- 
craft, by L. R. Woodworth and C. C. 
Kelber, Institute of the Aeronautical 
Sciences, Preprint No. 345, 1951. 

On the Theory of the Stationary and 
Pulsating Ramjet Engine (in German), 
by E. Sanger, Schweiz. Arch., vol. 16, 
November-December, 1950, pp. 11, 341- 
352, 369-378. 


Rocket Propulsion 
Engines 


Determination of Rocket-Motor Heat- 
Transfer Coefficients by the Transient 
Method, by Stanley Greenfield, Journal 
of the Aeronautical Sciences, vol. 18, 
August, 1951, p. 512. 

Heat Transfer in Rocket Motors and 
the Application of Film and Sweat Cooling, 
by R. H. Boden, Transactions of The 
American Society of Mechanical Engi- 
neers, vol. 73, May, 1951, pp. 385-390. 

Rocket Linings, by W. H. Duckworth, 
Aviation Week, vol. 54, June 25, 1951, 
p. 31. 

Survey of Development of Liquid 
Rockets in Germany and Their Future 
Prospects, by W. von Braun, Journal of the 
British Interplanetary Society, vol. 10, 
March, 1951, pp. 75-80. 

The Effect of Solid Particles in Rocket 
Exhausts, by N. J. Bowman, Journal of 
Space Flight, vol. 3, June, 1951, pp. 1-15. 

The Use of Auxiliary Rockets in High- 
Speed Aircraft, by Zb. Plaskowski, Air- 
craft Engineering, vol. 23, March, 1951, 
pp. 72-75. 


Heat Transfer and Fluid 
Flow 


Heat Transfer to a Fluid Flowing Tur- 
bulently in a Smooth Pipe With Walls at 
Constant Temperature, by R. A. Seban 
and T. T. Shimazaki, Transactions of The 
American Society of Mechanical Engineers, 
August, 1951, pp. 803-809. 

A Study of the Mechanism of Boiling 
Heat Transfer, by W. M. Rohsenow, 
Transactions of The American Society of 
Mechanical Engineers, vol. 73, July, 1951, 
p. 609. 

A Summary of Viscosity and Heat- 
Conduction Data for He, A, He, Oo, CO, 
CO,, H.O and Air, by F. G. Keyes, Trans- 
actions of The American Society of Mechani- 
cal Engineers, vol. 73, July, 1951, p. 597. 
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Investigation of Turbulent Flow and 
Heat Transfer in Smooth Tubes, Includ- 
ing the Effects of Variable Fluid Proper- 
ties, by R. G. Deissler, Transactions of 
The American Society of Mechanical Engi- 
neers, vol. 73, February, 1951, pp. 101- 
107. 

Fluid Flow through Porous Metals, by 
Leon Green, Jr., and Pol Duwez, Journal 
of Applied Mechanics, vol. 18, March, 
1951, pp. 39-45. 

Flow of Gases in Porous Media, by 
by L. H. Wilson, W. L. Sibbitt, and M. 
Jakob, Journal of Applied Physics, vol. 
22, August, 1951, pp. 1027-1030. 

Vaporization Rates and Heat Transfer 
Coefficients for Pure Liquid Drops, by 
Robert D. Ingebo, NACA _ Technical 
Note 2368, July, 1951. 

The Rate of Evaporation of Droplets. 
V. Evaporation Characteristics of Some 
Branched-Chain Hydrocarbons and a 
Straight-Chain Fluorocarbon, by R. §. 
Bradley and G. C. S. Waghorn, Proced- 
ings of the Royal Society of London, Series 
A, vol. 206, March 22, 1951, pp. 65-72. 

Some Kinematic Properties of Munroe 
Jets, by W. M. Evans, Research, vol. 3, 
August, 1950, pp. 376-378. 

Aerodynamic Heating at High Speeds, 
by T. Nonweiler, Journal of the British 
Interplanetary Society, vol. 10, July, 1951, 
pp. 160-176. 

Boundary Layers and Skin Friction in 
High-Speed Flow, by A. D. Young, 
Journal of the Royal Aeronautical Society, 
vol. 55, May, 1951, pp. 285-301. 

Investigations Concerning Flow Condi- 
tions in a Centrifugal Pump, and the 
Effect of Blade Loading on Head Slip, by 
J. F. Peck, The Institution of Mechanical 
Engineers, vol. 164, 1951, p. 1. 

Solution of Certain Unsteady Heat 
Flow Problems by Relaxation Methods, 
by A. Gilmour, British Journal of Applied 
Physics, vol. 2, July, 1951, pp. 199-204. 

An Air-Flow Analogy for the Solution 
of Transient Heat Conduction Problems, 
by M. B. Coyle, British Journal of Applied 
Physics, vol. 2, January, 1951, pp. 12-18. 

Compressible Flows, by J. D. Kelley, 
American Journal of Physics, vol. 19, 
April, 1951, pp. 237-243. 

Note on Desirability of Completely 
Expanded Nozzles, by C. B. Cohen and 
P. J. Evans, Jr., Journal of the Aeronau- 
tical Sciences, vol. 17, December, 1950, 
pp. 811-812. 


Combustion 


Net Heat of Combustion of AN-F-58 
Aircraft Fuels, by S. Rothberg and R. 8. 
Jessup, Industrial & Engineering Chem- 
istry, vol. 48, April, 1951, pp. 981-985. 

Method for Determining Distribution 
of Luminous Emitters in the Cone of a 
Laminar Bunsen Flame, by T. P. Clark, 
NACA Technical Note 2246, 1951, 28 
pp. 

Method of Calculating the Combustion 
Temperature of Gas Mixtures, by H. Wolff, 
Zeitschrift fiir Elektrochemie, vol. 54, 1950, 
pp. 556-560. 


Explosive Sensitivity of Ammonium 
Nitrate-Hydrocarbon Mixtures, by M. A. 
Cook and E. L. Talbot, Industrial and 
Engineering Chemistry, vol. 48, May, 1951, 
pp. 1098-1102. 

Thermodynamics and Chemical Kinet- 
ics of One-Dimensional Nonviscous Flow 
Through a Laval Nozzle, by 8. 8S. Penner, 
Journal of Chemical Physics, vol. 19, 
1951, pp. 877-881. 

Detonation of Secondary Explosives by 
an Electric Spark (in French), by R. 
Schall, Comptes rendus, vol. 233, 1951, pp. 
706-708. 

Elementary Processes Occurring in 
the Boudouard-Reaction (i.e., 2CO = 
CO. + Csotia + 41K cals) and Burning of 
Coal at Elevated Temperatures (in Ger- 
man), by H. Behrens, Naturwissenschaf- 
ten, vol. 8, 1951, p. 187. 

Chemiionization in Flames (Ionization 
induced by excessive concentrations of 
free radicals) (in German), by H. Beh- 
rens, Naturwissenschaften, vol. 8, 1951, 
pp. 187-188. 

On the Kinetics of Establishing of 
Equilibrium in Flames (in German), by 
H. Behrens, Zeitschrift ftir Physikalische 
Chemie, vol. 195, 1951, pp. 329-336. 

On the Structure and Combustion 
Kinetics of Diffusion Flames (in German), 
by H. Behrens, Zeitschrift fiir Physikal- 
ische Chemie, vol. 197, 1951, pp. 1-5. 

Self-Excitation of Flame Oscillations in 
Cylindrical Tubes (in German), by H. 
Behrens, Zeitschrift fiir Physikalische 
Chemie, vol. 197, 1951, pp. 6-16. 

Initiatium of Explosions in Liquids. 
Measurement of the Transient Pressures 
During Impact, by P. Gay, Transactions 
of the Faraday Society, vol. 46, pp. 848-852. 

The Theory of Flame Propagation I], 
by J. O. Hirschfelder and C. F. Curtiss, 
The Journal of Physical and Colloid Chem- 
istry, vol. 55, 1951, pp. 774-788. 

Physicochemical Aspects of the Com- 
bustion of Solid Fuels, by J. R. Arthur 
and J. R. Bowring, Industrial and Engi- 
neering Chemistry, vol. 43, 1951, pp. 528- 
533. 

Burning Velocities of Propane-Air Mix- 
tures, by F. H. Garner, R. Long, and G. 
K. Ashforth, Fuel, vol. 30, 1951, pp. 63-66. 

A Comparison of Quenching Distance 
and Inflammability Limit Data for Pro- 
pane-Air Flames, by Dorothy Martin 
Simon, Journal of Applied Physics, vol. 
22, 1951, p. 103. 

Pre-Combustion Chemical Reactions 
in Engines, by O. Widmaier, Brennstaf- 
Chemie, vol. 32, 1951, pp. 104-110. 

Ignition of Methane-Air Mixtures by 
Rapid Compression, by A. C. Edgerton, 
N. P. W. Moore, and W. T. Lyn, Nature, 
vol. 167, 1951, pp. 191-192. 

Study of the Flame-Stability Character- 
istics of Different Vaporized Fuels in a 
Small-Scale Combustion Chamber, by 
R. J. Wakelin, R. Heran, and R. R. Bald- 
win, Fuel, vol. 30, 1951, pp. 82-91. 

Flame Speeds in Hydrazine Vapour and 
in Mixtures of Hydrazine and Ammonia 
with Oxygen, by R. C. Murray, Trans- 
actions of the Faraday Society, vol. 47, 
July, 1951, pp. 743-751. 
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The Oxidation, Decomposition, Igni- 
tion, and Detonation of Fuel Vapors and 
Gases, by R. O. King, Canadian Journal 
of Technology, vol. 29, August, 1951, pp. 
382-390. 


A Discussion on Detonation, by W. G. 
Penney and others, Proceedings of the 
Royal Society of London, Series A, vol. 
203, November, 1950, pp. 1-33. 


The Cool Flame of Methyl Ethyl Ke- 
tone, by J. Bardwell and Sir Cyril Hinshel- 
wood, Proceedings of the Royal Society of 
London, Series A, vol. 205, February 22, 
1951, pp. 375-390. 


Spectroscopic Studies of Low-Pressure 
Flames V. Evidence for Abnormally High 
Klectronic Excitation, by A. G. Gaydon 
and H. G. Wolfhard, Proceedings of the 
Royal Society of London, Series S, vol. 
205, January 22, 1951, pp. 118-134. 


Flame Spectra of Systems with NO? 
Groups (in German), by H. Behrens and 
F. Rossler, Zeitschrift fiir Naturforschung, 
vol. 6a, 1951, pp. 154-160. 


Spectra and Kinetics of NO-Containng 
Flames (in German), by H. Behrens, 
Zeitschrift fiir Elektrochemie und Ange- 
wandte Physikalische Chemie, vol. 54, 
1950, pp. 535-538. 

The Emission of Flame Spectra, by A. 
G. Gaydon, Quarterly Reviews of the Chemi- 
cal Society of London, vol. 4, 1950, pp. 1-19. 

Kinetics of OH Radicals From Flame 
Emission Spectra III. Total Transition 


Probabilities and the Energetic Distribu- 
tion of OH(?2*) and O,(#2,,-) in the Oxy- 
Hydrogen Flame, by K. E. Shuler, Jour- 
nal of Chemical Physics, vol. 19, July, 
1951, pp. 888-894. 

Elementary Reaction in the Gas Phase 
Involving Excited Electronic States, by 
K. J. Laidler and K. E. Shuler, Chemical 
Reviews, vol. 48, April, 1951, pp. 153-224. 


Fuels, Propellants, and 
Materials 


Kinetics of the Polymorphic Changes in 
NH,NO; General Character of the Ki- 
netics of the IV-III Conversion, by B. V. 
Erofeev and N. I. Mitskevich, Zhur. 
Fiz. Khim., vol. 24, 1950, pp. 1235-1251. 

Stabilizer for Hydrogen Peroxide, by 
J. Isemura and S. Nakamura, Science 
(Japan), vol. 19, 1949, pp. 379-380. 

Force Constants, Frequency Assign- 
ments, and Thermodynamic Properties of 
Ozone, by M. J. Klein, F. T. Cleveland, 
and A. G. Meister, Journal of Chemical 
Physics, vol. 19, August, 1951, p. 1068. 

The Freezing-Point Diagram of the 
Hydrazine-Water System, by T. C. H. 
Hill and J. F. Sumner, Journal of Chemical 
Society, March, 1951, pp. 838-40. 

The Heat of Dissociation of Fluorine, 
by R. N. Doescher, Journa! of Chemical 
Physics, vol. 19, August, 1951, p. 1070. 


Autoxidation of Hydrazine, by L. F. 
Audrieth and P. H. Mohr, Industrial and 
Engineering Chemistry, vol. 43, August, 
1951, p. 1774. 

Ignition Temperatures of Magnesium 
and Magnesium Alloys, by W. M. Fassell, 
Journal of Metals, vol. 3, July, 1951, p. 
522. 

The Manufacture of Sintered Alumina 
Ware, by R. Scott, Metallurgia, vol. 43, 
July, 1951, pp. 50-52. 

The Preparation of Alumina Bodies by 
the Slip-Casting Method, by W. Watt and 
J. P. Roberts, Metallurgia, vol. 43, June, 
1951, pp. 307-308. 


Engineering Materials in Ordnance 
Equipment, by B. S. Mesick, Materials 
and Methods, vol. 33, May, 1951, p. 59. 

Titanium and Its Alloys, by Anderson 
Ashburn, American Machinist, vol. 95, 
June 11, 1951, p. 146. 

Linear Thermal Expansion of Arti- 
ficial Graphites to 1,370 C., by M. D. 
Burdick, Journal of Research, vol. 47, 
July, 1951, pp. 35-40. 


Some Mechanical Properties of Graphite 
at Elevated Temperatures, by C. Malm- 
strom, R. Keen, and L. Green, Jr., Journal 
of Applied Physics, vol. 22, May, 1951, 
pp. 593-600. 

An Introduction to Arc-Cast Molybde- 
num and Its Alloys, by J. L. Ham, Trans- 
actions of The American Society of Me- 
chanical Engineers, vol. 73, August, 1951, 
pp. 723-732. 


Ignition Analyzers, 


SCINTILLA MAGNETO DIVISION 


BENDIX AVIATION CORPORATION 


SIDNEY, NEW YORK 


Manufacturers of Ignition Systems for Jet, Turbine and 
Piston Power Plants, Fuel Injection Equipment for Railway, 


Marine and Industrial Diesel Engines, Electrical Connectors, 


Moldings and Special 


Components. 


NovemMBER, 1951 


195 


* 
md 
51, 
1et- 
low 
ner, 
19, 
by 
PP. 
in 
in 
haf- 
tion 
of a 
Beh- 
951, 
of 
ische 
tion 
1an), 
ikal- 
ns in a 
H. 
ische 
uids. 
tions 
-852. 
“4 

acter 
in a 
Bald- 


High Temperature Steels and Alloys for 
Gas Turbines, by T. Bishop, Metal Prog- 
ress, vol. 59, May, 1951, pp. 653-656, 
690-696. 

Some Properties of Beryllium Oxide 
and Beryllium Oxide-Niobium Ceramals, 
by C. F. Robards and J. J. Gangler, 
NACA Research Memo E50G21, March, 
1951, 18 pp. 


Physical-Chemical Topics 


Availability and Irreversibility in Ther- 
modynamics, by Joseph H. Keenan, 
British Journal of Applied Physics, vol. 2, 
July, 1951, pp. 183-192. 

Irreversibility and Generalized Noise, 
by H. B. Calhn and T. A. Welton, Physicals 
Review, vol. 83, 1951, pp. 34-40. 

Viscosity Behavior of Gases, by L. A. 
Bromley and C. R. Wilke, Industrial and 
Engineering Chemistry, vol. 43, July, 
1951, pp. 1641-1648. 

Symposium on Anomalies in Reaction 
Kinetics—Four Sessions, Journal of Phys- 
tcal and Colloid Chemistry, vol. 55, 
June, 1951. 

Air Condensation in Hypersonic Flow, 
by P. Wegener, S. Reed, Jr., E. Stollen- 
werk, and G. Lundquist, Journal of Ap- 
plied Physics, vol. 22, August, 1951, pp. 
1077-1085. 

Recent Theories of Transport Processes 
in Liquids, by R. Eisenschitz, Nature, 
vol. 167, February, 1951, pp. 216-220. 


Instrumentation and 
Experimental Techniques 


A Method for Correcting Orifice-Meter 
Measurements for Flow Pulsations, by E. 
V. Hardway, Jr., Instruments, vol. 24, 
July, 1951, pp. 763-765. 

A New Flash Illumination Unit for 
Ballistic Photography, by H. F. Quinn 
and O. J. Bourque, Review of Scientific 
Instruments, vol. 22, February, 1951, pp. 
101-105. 

A Study of Head Loss in Venturi-Meter 
Diffuser Sections, by Joel Warren, T'rans- 
actions of The American Society of Mechani- 
cal Engineers, vol. 73, May, 1951, pp. 399- 
402. 

An Analog Computer for Flame Gas 
Composition, by W. S. McEwan and §. 
Skolnik, Review of Scientific Instruments, 
vol. 22, March 1951, pp. 125-133. 

An Electromagnetic Flowmeter for 
Transient Flow Studies, by J. S. Arnold, 
Review of Scientific Instruments, vol. 22, 
January, 1951, pp. 43-47. 

Discharge Measurements by Means of 
Venturi Tubes, by A. L. Jorissen, Trans- 
actions of The American Society of Me- 
chanical Engineers, vol. 73, May, 1951, pp. 
403-411. 

Measurement of Velocity and Pressure 
of Gases in Rocket Flames by Spectro- 
scopic Methods, by F. P. Bundy, H. M. 
Strong, and A. B. Gregg, Journa! of Ap- 
plied Physics, vol. 22, August, 1951, pp. 
1069-1077. 

High-Speed Sampling Techniques, by 
B. R. Shepard, Electronics, vol. 24, August, 
1951, pp. 112-115. 
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proving the routine practicality of flying far beyond the 
speed of sound. As in the Air Force’s historic X-l, the 
Skyrocket’s successful rocket engine was designed and built 
by REACTION MOTORS. 


The broadest background of research and development in 
the field have produced an incomparable record of rock- 
etry achievements for Reaction Motors, Inc. ... progress so 
continuous that yesterday’s rocket engines already are 
obsolete. If you haven‘t talked with RMI today, you do not 
know the newest developments in rocket power for piloted 
aircraft, guided missiles and many other applications. 


~ 


SINCE 1941—THE FIRST NAME IN RQCKETRY 
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field of the future. Design and layout draftsmen—stress analysts— 
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Induction Method for Measuring Fluid 
Flow, by A. J. Morris, Electrical Engi- 
neering, vol. 70, June, 1951, p. 529. 

Instruments for Upper Atmosphere and 
Interplanetary Navigation, by J. C. 
Bellamy, Navigation, vol. 2, December, 
1950, pp. 272-275. 

Kerr Cell Photography of High Speed 
Phenomena, by E. M. Pugh, R. v. Heine- 
Geldern, S. Foner, and E. C. Mutschler, 
Journal of Applied Physics, vol. 22, April, 
1951, pp. 487-493. 

Power Unit for High-Intensity Light 
Source, by A. E. Young, 8. McCullough, 
and R. L. Smith, NACA Research Memo 
150K27, July 9, 1951. 

Solar Furnace for High Temperature 
tesearch, by W. M. Conn, American Ce- 
ramic Society Bulletin, vol. 30, 1951, pp. 
163-164. 

Step Multiplier in Guided Missile Com- 
puter, by Edwin A. Goldberg, Electronics, 
vol. 24, August, 1951, pp. 120-124. 

Study of Linear-Resistance Flowmeters, 
by F. W. Fleming, 7’ransactions of The 
American Society of Mechanical Engi- 
neers, vol. 73, July, 1951, p. 621. 

The Long Range Proving Ground, by 
Robert D. Hullihan, Ordnance, vol. 36 
July-August, 1951, pp. 158-160. 

The Trajectograph ONERA (in French), 
by I. Podiasky, La Recherche Aéronau- 
tique, Nov.—Dec., 1950, pp. 35-48. 

Use of Image Converter Tube for High- 
Speed Shutter Action, by A. W. Hogan, 


Proceedings of the Institute of Radio Engi- 
neers, vol. 39, March, 1951, pp. 268-270. 


Woomera; Post-War Progress with 
British Guided Missiles; The Australian 
Rocket Range, Flight, vol. 59, April 13, 
1951, pp. 429-431. 


Flight, Ballistics, and 
Vehicle Development 


Research with High-Speed (Rocket- 


Powered) Models, Anonymous. Aero- 
plane, vol. 81, August 17, 1951, pp. 190- 
192. 


Fundamental Dynamics of Reaction- 
Powered Space Vehicles, by L. N. Thomp- 
son, Aircraft Engineer, vol. 23, August, 


1951, p. 228. 


The Evasion of Hazardous Objects in 
Space, by Wayne Prvell, Journal of Space 
Flight, vol. 3, May, 1951, pp. 1-9. 

Meteors and Space Travel, by Michael 
W. Ovenden, Journal of the British Inter- 
planetary Society, vol. 10, July, 1951, 
pp. 176-180. 


Evolution of the Guided Missile, by 
K. W. Gatland, (Series beginning May 4th) 
Flight, vol. 60, August 3, 1951, pp. 140-143 

Motion Relative to the Surface of the 
Rotating Earth, by R. H. Bacon, Ameri- 
can Journal of Physics, vol. 19, January, 
1951, nn. 52-56. 


Shock and Vibration 
Measurements 


Mechanical 
Electrical Research 


of 


Telemetering 
Applications 


Geophysical 
Explorations 
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On the Accuracy of Long-Range Ballis- 
tic Rockets, by W. E. Frye, Journal of 
Applied Physics, vol. 22, May, 1951, pp. 
585-589. 

What Will Space Flight Require? 
Eugen Sanger, Rocketscience, vol. 5, June, 
1951, pp. 26-33. 

The Lunar Base, by G. V. E. Thomp- 
son, Journal of the British Interplanetary 
Society, vol. 10, March, 1951, pp. 49-70. 

Long-Distance Rocket Flight, Anony- 
mous, Mechanical Engineering, vol. 73, 
August, 1951, p. 656. 

Orbital Rockets: I Some Prelimi- 
nary Considerations, by Kenneth W. Gat- 
land; II The Rocket Structure with Spe- 
cial Reference to Expendable Construc- 
tion, by Alan E. Dixon; III Conception of 
an Instrument-Carrying Orbital Rockei, 
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Jet Engine 
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by Anthony M. Kunesch, Journal of the 
British Interplanetary Society, vol. 10, 
May, 1951, pp. 97-115. 


Other Classifications 


(Meteorology, Astrophysics, etc.) 


A Preliminary Study of a Physical Basis 
of Bird Navigation, by H. L. Yeagley, 
Journal of Applied Physics, vol. 22, June, 
1951, pp. 746-760. 


The Effects of Time and Temperature 
Upon Tolerance to Positive Acceleration, 
by Erneste Martin and James P. Henry, 
Technical Data Digest, vol. 16, June, 1951, 
pp. 19-23. 


Physics and Psychophysics of Weight- 
lessness, by H. Haber, Journal of Aviation 
Medicine, June, 1951, p. 180. 


Influence of Certain Anti-Motion-Sick- 
ness Drugs on Psychomotor and Mental 
Performance, by R. A. MacKay, Journal 
of Aviation Medicine, June, 1951, p. 194. 


Biological Problems of Space Flight, by 
J. B. Haldane, Journal of the British Inter- 
planetary Society, vol. 10, July, 1951, pp. 
154-158. 


Physiological Endurance in Jet Missiles 
(in Italian), by G. A. Crocco, Termotec- 
nica, vol. 5, February, 1951, pp. 55-59. 


Radiation Effects on Man in Space, by 
Conrad Buettner, American Meteor So- 
ciety Bulletiin, May, 1951, p. 183. 


ROCKETS 


by 


Dr. R. H. Goddard 


A theoretical and experimental classic by 
the rocket pioneer. Includes “A Method of 
Reaching Extreme Altitudes’ and “Liquid 
Rocket Propellant Development.” 1946, 111 
pages. Price $3.50. Order from the American 
Rocket Society, 29 West 39th Street, R-719, 
New York 18, N. Y. 
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The Air Force Missile Test Center, used by all our 
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over the Bahamas, into the South Atlantic. 
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air supremacy! 
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Air Force B-61 Matad 
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- the first to use the completed range. THE Gienn L. 
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of research and development of jet and rocket propulsion devices 
and their application to problems of transportation and com- 
munication. At the present time it is also concerned with the 
military uses of the reaction principle. 
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Aerojet’s liquid 


and solid assist-takeoft 


rockets and new, improved, 


forward firing aircraft 
rockets are an important 
part of America’s 


air defense. 


Liquid and solid rockets 
for assist-takeoff. 
Sounding rockets and 
missile boosters. 
Rocket test 
installations. 

Complete research, 
engineering, and 
manufacturing 
facilities. 
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